


7 


` QA 


; N N 
kg , ! ts | 









Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1975 


A study of unsteady transonic interference effects. 


Schlein, Paul Barron 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/21069 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
P (q D U DLEY research materials and institutional publications created by the NPS community. 
FW | 3 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
" ll KNOX appointed — and published — scholarty author. 
http://www.nps.edu/library 






LIBRARY 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 











- Gum | ; Some > 


ATA PE DOI a BPE LAO 








INTERFERENCE EFFECT 
by 
Paul Barron Schleirn 


March 1975 














DESTIDY OF UNSTEADY TRANSONIC 
io o 





A vos. OTRA 


Thesis Advisor: MIT Patzer 
v 242 EEE TAT ANETTE TAU EE, RETTEN ITO LIT ἽΝ IL 
Approved for public release; distribution untimited. 





UNCLASSIFIED 


SS “πατα επι -πΏ-εππακώκνταπα.οἹ---ἠρ-πορ οφ --ν-----ο---ν για οπιαικκκκιασσημώνετ ίσα επικ, ιαερικαρτρκορ (μακαρι. 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE 


. REPORT NUMBER 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


2. GOVT ACCESSION NO. 3. RECIPIENT'S CATALOG NUMBER 









S. TYPE OF REPORT & PERIOD COVERED 


Doctor of Philosophy 


March 197 
6. PERFORMING ORG. REPORT NUMBER 


4. TITLE (and Subtttle) 





A Study of Unsteady Transonic 
Interference Effects 











7. AUTHOR(e) 


Paul Barron Schlein 


8. CONTRACT OR GRANT NUMBER(e) 









9. PERFORMING ORGANIZATION NAME ANO ADORESS 10. PROGRAM ELEMENT, PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


Naval Postgraduate School 
Monterey, California 93940 


Σι. CONTROLLING OFFICE NAME AND ADDRESS 12. REPORT DATE 
March 1 
Naval Postgraduate School 75 μον. 
Monterey, California 93940 176 


4. MONITORING AGENCY NAME & ADORESS(iIf different from Controlling Oflice) 15. SECURITY CLASS. (of thie report) 





Unclassified 


15a, OECL ASSIFICATION/ DOWNGRADING 
SCHEOULE 


16. DISTRIBUTION STATEMENT (of thle Report) 


Approved for public release; distribution unlimited. 


17. OISTRIBUTION STATEM ENT (of the abstract entered In Block 20, ff different from Report) 


18. SUPPLIMENTARY NOTES 


| 18. ΚΕΥ WOROS (Couvltnue on reverso etde If neceeeery and identify by dfock number) 
Unsteady Aerodynamics 
Aeroelasticity 
Cascade Aerodynamics 
Transonic Aerodynamics 
Wind Tunnel Interference 
[20. ABSTRACT (Gontinue on reveres efde ff neceesery end identity by block number) 

Tne unsteady, sonic, small disturbance potential equation was approximated 
by reduction to parabolic form and sclved by Laplace and Fourier Transform and 
Collocation techniques, including thickness effects through the Oswatitsch 
parameter. The Transform techniques produced reliable solutions and the 
Collocation solution was usable with qualifications, for some parameter values, 
stemming from the choice of approximating series. Stability boundaries for an 
unstaggered cascade oscillating in sonic flow were calculated with the Fourier 


———"————— A A O DE A E 


DD - δῶν 2 1473 EDITION OF ! NOV $5 IS OBSOLETE UNCLASSIFIED 
(Page 1) Š | SECURITY CLASSIFICATION OF THIS PAGE (When Date En'eret) 








UNCLASSIFIED 


a A eng u | u — = 
SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 


(20. ABSTRACT Continued) 





Transform technique revealing that blade interference is strongly destabiliz- 
ing. The Collocation technique was successfully applied to the two- 
dimensional, unsteady, porous walled, sonic wind tunnel interference problem, 
heretofore unapproached with the inclusion of thickness effects. 





DD Form 1473 @ACK) > 
LS i UNCLASSIFIED 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(When Dats Entered) 





A Study of Unsteady Transonic 
interference Effects 


by 


Paul Barron Schlern 
Lieutenant, United States Navy 
B.S., United States Naval Academy, 1966 


Submitted in partial fulfiilment of the 


requirements for the degree of 


DOCTOR OF EHILOSOPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1975 





ABSTRACT 


The unsteady, sonic, small disturbance potential equation 
was approximated by reduction to parabolic form and solved by 
Laplace and Fourier Transform and Collocation techniques, 
including thickness effects through the Oswatitsch parameter. 
The Transform techniques produced reliable solutions and the 
Collocation solution was usable with qualifications, for some 
parameter values, stemming from the choice of approximating 
series. Stability boundaries for an unstagpered cascade 
oscillating in sonic flow were calculated with the Fourier 
lp sermotecehnique revealing that blade interference 15 
EHmsonelyodestabilrzing. The Collocation technique was suc- 
cessfully applied to the two-dimensional, unsteady, porous 
walled, sonic wind tunnel interference problem, heretofore 


unapproached with the inclusion of thickness effects. 





Er. 


DIL; 


TABLE OF CONTENTS 


πα O A a 15 
UNSTEADY TRANSONIC FLOW THEORY -~--------------------- 21 
A. EQUATIONS FOR GENERAL OSCILLATORY MOTION -------- 24 
B. REDUCTION OF THE GOVERNING EQUATION ------------- 22 
που: οσο”τιοῖ”.....------------------------- An 

1. Flow Tangenoey ------------------------------- 2l 

2. Zero Upstream Disturbance ------------------- 25 
FORMULATION OF SOLUTIONS ---------------------------- 26 
A. SOLUTION THROUGH THE LAPLACE TRANSFORM ----------- 26 

1. Cascade Geometry and Boundary Conditions ---- 26 


ου. ουσ οαῃ of the Oscillatory Sonic 


Flow Equation ------------------------------- 28 

3. Modification for Sonic Free Jet Problem ----- 32 

B. SOLUTION THROUGH THE FOURIER TRANSFORM ---------- 33 
1. Cascade Geometry and Boundary Conditions ---- 33 
ans iormabtion of the Sonic Flow Equation --- 35 

3. Inversion ----------------------------------- 38 


h. Calculation of Lift and Moment Coefficients - HO 


eee OOlutions for Sy and Tw ———————-——----------- 41 
C. SOLUTION THROUGH COLLOCATION TECHNIQUE ---------- 45 
D. EVALUATION OF THE INTEGRALS OCCURRING IN 

EQUATIONS (III-24,86,87 and 88)------------------- 51 

l. The Steady Case ----------------------------- 52 

2. The Oscillatory Case with A = 0 ------------- 5 





Hy 


Vi. 


3. The General Case, K Z O and À Z 0 ----------- 56 
4, The Case for P = 0, Jj = 2, K #0 and A #0 -- 56 


E. APPLICATION AND COMPARISON OF SOLUTIONS TO 


SEVERAL CASCADE PROBLEMS ------------------------ 56 
l.  Interblade Phase Angle Variation ------------ 50 
2. Frequency and ος Variation ----------- 59 
3. Oswatitsch Parameter = 0.0 ------------------ 70 
TRANSONIC CASCADE AND BIPLANE STABILITY ------------- 76 
A. GENERAL ----------------------------2------------- 76 
B. TRANSONIC, UNSTAGGERED CASCADE STABILITY -------- 77 
C. TRANSONIC BIPLANE STABILITY --------------------- 83 
ENTENSTON” OE THE COLLOCATION SOLUTION TO THE 
POROUS WALLED TRANSONIC WIND TUNNEL ----------------- 85 
A. PROBLEM FORMULATION ----------------------------- 85 
l. Geometry and Boundary Conditions ------------ 85 


PERSO TuUt ton of the Oscillatory Sonic Flow 
Equation ----------------------------------------------------------------------- 85 


3. Derivation of Expressions for Derivatives 
of the Potential ---------------------------- 86 


B. RESULTS OF POROUS WALL WIND TUNNEL CALCULATIONS - 90 


C. DISCUSSION OF MURASAKI AND DRAKE POROUS WALL 


SONIC WIND TUNNEL SOLUTIONS --------------------- 08 
DISCUSSION AND SUMMARY ------------------------------ 101 
A. ANALYSIS OF SOLUTIONS --------------------------- 101 

1. Laplace Transform m. --------------------- τοι 

2. Fourier Transform --------------—------------- ΠΠ 

3: Collocation Νετποᾶ -------------------------------.-γ----- T02 
B. CASCADE STABILITY ------------------------------- 103 





C.  POROUS WALL WIND TUNNEL INTERFERENCE ------------ 101 
D. FUTURE WORK ------------------------------------- 104 
VII. CONCLUSIONS -------------------------------2---------- 207 
APPENDIX A: LAPLACE TRANSFORM SOLUTION COMPUTER ROUTINE --109 


APPENDIX B: FOURIER TRANSFORM SOLUTION COMPUTER ROUTINE --126 


APPENDIX C: COLLOCATION TECHNIQUE COMPUTER ROUTINE ------- 133 
APPENDIX D: POROUS WALL WIND TUNNEL COMPUTER ROUTINE ----- 153 
E O AA 8. l l... l 2... =... 172 


INITIAL DISTRIBUTION LIST ------------ ------ - ------ - -- ----- I9 


-- 





LIST OF TABLES AND FIGURES 


TABLES 


n. COMPARISON OF POTENTIAL (REAL/IMAGINARY) CALCU- 
LATIONS OF COLLOCATION, FOURIER, AND LAPLACE 
SOLUTIONS FOR VARIOUS INTERBLADE PHASE ANGLE ------- 


MEGURES 
d. UNS@AGGE RED CASCADE ---2------------------------.----- 
p. UMNSBNCGEREDOCASCADES ----2---------------------------- 


ga-T. COMPLEX VELOCITY POTENTIAL COMPARISONS FOR o = T 
AND p = 1.0 WITH VARYING FREQUENCY AND OSWATITSCH 
PARAMETER. VELOCITY POTENTIAL VS X ---------------- 


la-d. COMPLEX VELOCiITY POTENTIAL COMPARISONS FOR oc - T 
AND À = O, VARYING FREQUENCY AND INTERBLADE 
BER DEN ICE VELOCITY POTENTIAL VS Ὁ ----------------- 


Sa-e. CASCADE STABILITY BOUNDARY FREQUENCY VS PIVOT 
E ο στ ο JJ ποστς 
5-f. BIPLANE STABILITY BOUNDARY FREQUENCIES VS PIVOT 
EIE NEL A Z. e e e ee LLL 
ο, ΤΕ ΡΕ PLOW PRESSURE COEFFICIENT VS X, t = +06, 
πο eee hOUS POROSITY (0. ,.05,.50,10. ,100.) AND 
SOLID WALL TUNNEL, SONIC FREE JET, AND EXPERIMENTAL 
RESULTS -------------------------------------------- 


7a-b. NONSTEADY FLOW NORMALIZED PRESSURE COEFFICIENT 
Der = .2, 1 = .06, FOR VARIOUS POROSITY 
(0.,.05,.50,5.0,100.), SOLID WALL TUNNEL, SONIC 
FREE JET, AND VARIABLE POROSITY COLLOCATION 
ROUTEN Husa a ee ee = 


8a~b. NONSTEADY FLOW NORMALIZED PRESSURE COEFFICIENT 
VS X, K = .5, T = .06, FOR VARIOUS POROSITY 
(0. ,.05,1.0.5.0,100.), SOLID WALL TUNNEL AND 


ο ο ρα JET RESULTS ------2- _- -_ _ ZZ 222... __ 





Exp[ a] 
E 


E; 


LIST OF SYMBOLS 
K^/(A* + uk?) 


doned mensional position of pivot location, 
a divided by reference length 


postition or pivot in-pitch oscillation 


H orlan lemor ler Transform solution, functions 
of interblade phase angle or can be = AP*/Y 


nondimensional pivot location, b divided by 
reference length 


airfoil or blade semichord 


dimensionless speed of' sound; speed of sound 
divided by reference velocity 


ji e meo í(r'olent y positive upwara 

moment coefficient, positive leading edge up 
pressure coefficient 

wn Oil or blade chora 


variable in Fourier Transform solution, function 
of interblade phase angle 


a 
€ 


KP^/l 


constant in solution of Fourier transformed 
potential equation 


constant in solution of Fourier transformed 
potential equation 


e o entegral of Y, used in Fourier Transiorm 
solution 


η integrals or «5.1 7 4 


variable necessary in Fourier Transform 
solution = 1,2 





H(X,T) 


H(Z) 


5| 


n(X) 


I(X) 


J(X) 


menmdmenmsiosal Y position of airfoil or blade 
surface 


Heaviside step function 
nondimensional interblade distance, P/C 


nondimensional mean Y position of airfoil or 
blade surface 


nondimensional ampditude of oscillatory airfoil 
or blade surface . 


amplitude of plunge oscillation 


nondimensional amplitude of plunge oscillations, 
divided by reference length, l = 1 or 2 


integrals appearing in Laplace Transform and 
Collocation solution 


integral values for Fourier Transform solution 
imaginary constant, y-1 , or variable index 
integral values for Fourier Transform solution 
reduced frequency; ub/U. 

uc/U4 


Oswatitsch parameter, A, or number of collocation 
points to be used 


Mach number of reference velocity 
variable index 

interblade distance, or variable in La 
Laplace Transform transformation variable 
À + 21K 

real part 


variable of integration for Laplace Transform and 
Collocation solutions 


see Ty 


10 





dimensionless time, t multiplied by reference 
velocity and divided by reference length 


integral values necessary in Fourier Transform 
solution 


time 
X-S, variabie used in evaluation of I 
reference velocity, freestream or critical 


ven terence source distribution, m = 0 or 1 for 
reference or adjacent blade, respectively 


nondimensional downwash velocity at blade chordline, 
used in Laplace Transform and Collocation solutions 


source distribution on airfoil or blade, m = 0,1, 
for reference or adjacent blade, respectively 


nondimensional downwash velocity at blade chordline, 
Weeden Outer Ivansform solution 


dimensionless coordinate in freestream direction; 
x divided by reference length 


Ronrdanmensional pivot location, e divided by Ὁ 
AL 

upper limit of 'ceorrection integral! 
coordinate in freestream direction 


aimensionless coordinate perpendicular to stream 
direction 


coordinate perpendicular to stream direction 
amplitude of pitch oseillation 
ratio of specific heats 


perpendicular coordinate from blades, i=0,1 





variable in summation for J(X) 

coefficlents of isolated blade source distribution 
series, m defines blade (reference or adjacent), 

n defines sequence number in series 

Oswatitsch parameter, (y+1) dy, 

variable in transformed potential equations 
Fourier Transform transformation variable 


variable in summation for I(X) 


coefficients of interference source distribution 
series, see ϐ 
mn 


coordinate along airfoil or blade chord from 
midchord or leading edge, i = 1 or 2, respectively 


nondimensional interblade distance, P divided by b 
interblade phase angle 

porosity parameter 

airfoil or blade thickness/chord ratio 
Gimensionless unsteady two-dimensional velocity 
τυ}. potential divided by reference length 


and reference velocity 


nondimensional perturbation potential, or dummy 
variable 


isolated blade potential for the reference blade 
isolated blade potential for the adjacent blade 
interference potential for the reference blade 
interference potential for the adjacent blade 


Steady fiow perturbation potential, nondimensional 


T2 





«| 


nonsteady flow perturbation potential, nondimensional 
transformed nonsteady perturbation potential 


circular frequency of pitch or plunge oscillation 


T) 





ACKNOWLEDGEMENT 


This work could not have been completed were it not for 


the indefatigability of Dr. Max F. Platzer, the timely 


assistance of Dr. Craig Comstock, and the inspiration of 


Louisa May Alcott: 


Far 
are 
EID 
and 
and 


away in the sunshine 

my highest aspirations. I may 
reach them, but 1 can look up 

see their beauty, believe in them 
try to follow where they lead. 


11 





I. INTRODUCTION 


Nonsteady, interfering, transonic flows are of increasing 
importance to aeroelasticians. Recent design concepts for 
|psnsonie transports, high maneuverability fighters, biplanar 
Space shuttle configurations, and so called transonic axial 
flow compressors, to name a few, involve such flows when aero- 
elastic problems are considered.  Ironiealiy, while sub/super- 
sonic flutter problems have been extensively studied, the 
transonic problem has been circumvented by either rapid 
acceleration through the regime or operation avoiding it. It 
is now apparent that the question of the effect of transonic 
flow, including possible adjacent flow field interference, 
on flutter must be confronted in a practical way. 

Of course, even without interference, nonsteady transonic 
flows challenge the &eroelastician concerned with aerodynamic 
apie or structural vibrations. In that regime, pressure 
perturbations generated by a body and propagated at the local 
Speed of sound tend to accumulate near the body itself. Thus, 
the magnitude of aerodynamic forces is increased and, perhaps 
more importantly, a large phase lag can exist between the 
forces and the motion producing them. It is this phase lag 
that can lead to negative aerodynamic damping and, potentially, 
dangerous flutter. 

Πο  usenough to analytically study the flutter charac- 


πι Darvicular configurations. Wind tunnel experiments 


ies 





must be conducted to verify analytical results before a 
configuration ean be committed to flight test. However, 

flow perturbations from a wind tunnel test model may reflect 
trom the tunnel walls back onto the model requiring a correc- 
tion to test results due to interference. In steady transonic 
flow, such wall interference has been thoroughly examined and 
porous opr slotted wall tunnels have been utilized for some 
time to reduce its effect. In nonsteady flow, the porous 
uu boundary condition is a difficult one for classical 
analytical approaches to handle. Only Drake [Ref. 6] has in- 
cluded wall porosity in a theoretical analysis of a thin air- 
foil oscillating in a transonic wind tunnel and his results 
are by no means complete. Thus, there is a great need for 
solutions that can yield the magnitude of unsteady wall 
interference and the effect of wall porosity upon it. 

Blade flutter in transonic axial compressors is a complex 
problem involving, among other things, interference of highly 
perturbed transonic flow fields. Historically, sub and super- 
sonic axial compressor flows have been simulated by a series 
of airfoils (or blades) placed in close proximity (a cascade) 
so that the two dimensional flow models that in a rotating 
compressor stage. For the transonic compressor stage the 
mixed flow character of the real flow makes the two-dimen- 
sional approximation questionable, e.g., McCune [Ref. 16]. 

It must be recognized that accurate modeling of such a complex 


flow field is well beyond the state of the art. However, the 





cascade approximation may yield important information relative 
centhe trends of blade flutter. 

Nonsteady transonic flow analysis begins with the deriva- 
tion of the nonlinear small perturbation potential equation 
which Lin, Reissner, and Tsien [Ref. 15] and Landahl [Ref. 

13] have shown can be linearized for sufficiently high fre- 
quency. Since aeroelastic analyses are, generally, concerned 
with the stability of small amplitude oscillations, the solu- 
tions to the nonlinear potential equations shall be written 
as the sum of a steady component and a much smaller nonsteady 
component. The result is two equations, the usual nonlinear 
steady flow equation and a linear nonsteady flow equation 
containing variable coefficients depending upon the steady 
flow. 

Transonic cascade flow for thin blades and in-phase oscil- 
lations was approached by Hamamoto [Ref. 10] using Fourier 
Transform techniques on the linearized potential equation. 
Gorelov [Ref. 8] solved the linear nonsteady equation for 
mildly supersonic flow past thin blades with a superposition 
of isolated blade and interference potentials represented as 
trigonometric series. See also Samoylovich [Ref. 221. 

Elder [Ref. 7] revised Gorelov's Collocation technique by 
using a power series and solved for the steady sonic flow 
cascade potential and the coefficients of the series for the 
nonsteady flow, out-of-phase oscillation case. Elder made use 


of the Oswatitsch [Ref. 20] parameter (used by Oswatitsch for 





steady flows), as ka y Teipel [Ref. 30] and Hosokawa 
[Ref. 11] for nonsteady flows, to allow for blades of nonzero 
thickness (nonuniform steady potential flow). 

Sonic wind tunnel interference in steady flow was solved 
by Sandeman [Ref. 23] using Spreiter's local linearization 
technique, this after Guderley [Ref. 9] and Spreiter et al 
[Ref. 27] showed that the choked tunnel results can be cor- 
rected so as to simulate unbounded flows. Murasaki [Ref. 18] 
selved steady transonic tunnel fiows for porous or teca 
walls. Miles [Ref. 17] and Drake [Ref. 5] solved the super- 
sonic nonsteady solid and porous wall cases, respectively. 
Drake [Ref. 6], then, utilized Laplace Transform techniques 
and the linear potential equation for a thin airfoil between 
porous walls for transonic flows. Recently, Savkar [Ref. 24] 
Mas Laplaced transformed the linearized equation and studied 
roll stability in a solid walled wind tunnel. 

Wind tunnel interference can be approached, according to 
Sears [Ref. 25], with suitably located sensors within the 
test section (not on the model) and variable wall shape or 
porosity and/or plenum pressure behind the wall. The sensors 
would modify the tunnel parameters until unbounded flow con- 
ditions were present at them. Sears argued that unbounded 
flow conditions are thus obtainable upon the model. 

This work began with the object of applying a Collocation 
technique to the porous wall wind tunnel problem with the con- 


viction that that method could suitably handle the classically 


18 





' !τ1πιὶὐ porous wall boundary condition. While proceeding 
toward that goal, an efficient approximate solution was de- 
veloped for studying the aerodynamic stability of the sonic 
oscillating cascade and stability boundaries were calculated 
for the cascade as well as for a biplanar configuration. 

To accomplish the initial objective of this work, the 
nonsteady potential equation was reduced with the aforemen- 
tioned Oswatitsch parameter to approximate the variable 
coefficient therein. The EE reduced equation was 
solved by a Collocation technique for the cascade problem 
following Elder [Ref. 7]. For numerical checks upon the 
Collocation technique, the equation was Laplace transformed 
extending the steady flow work of Sandeman [Ref. 23] and 
Fourier transformed in a manner analogous to Hamamato [Ref. 
10]. These extensions were made flexible enough to handle 
blade thickness variations (through the Oswatitsch parameter) 
and variable interblade phase angle. 

The Elder power Series initially utilized in the Collo- 
cation solutions to the cascade problem was found to have 
numerical difficulties that increased computation time re- 
quirements and reduced accuracy. Thus, a 'step-wise-linear' 
series was introduced and the resulting Collocation method 
solution to the cascade problem agreed well enough with the 
r two solutions to allow Some confidence in its subsequent 


adaptation to the porous wall wind tunnel problem. 
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The Collocation solution was found to be quite flexible. 
It can be altered to allow streamwise variation of the Os- 
watitsch parameter, in a manner reminiscent of Spreiter's 
local linearization, to improve accuracy. Further, it allows © 
a variation in wall porosity for possible simulation of the 


technique mentioned by Sears. 
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INS TEADY TRANSONIC FLOW THEORY 


A. EQUATIONS FOR GENERAL OSCILLATORY MOTION 

The governing equation considered below arises from the 
assumption that the flows of interest are of frictionless 
perfect gas. Thence, the equation Of state and Eulerian 
equations of motion are invoked and simplified, assuming the 
absence of heat conduction and strong shock processes. Fur- 
ther, a velocity potential is introduced which reduces the 
number of dependent variables at the imposition of the re- 
quirement of flow irrotationality. Finally, a small distur- 
bance analysis, including the unique properties of transonic 


flows (see Refs. 12 and 13), yields, in two dimensions 
0 a (yO WG = eo, = 0. (TIN) 
Tr XT Mé p KX uc Dos j 


Mats equation iS applicable, approximately, to flows containing 
weak shocks, incidentally (see Ref. 13). 
In the study of unsteady aerodynamics, it is customary 
and convenient to expand the solution into steady and 
oscillatory nonsteady components. Thus, 
COG Vey a=) GY) ot R. P. [o (x, Y) e 17] (I1-2] 


3 


where K is the reduced frequency. 





Substitution in equation (II-1) yields a coupled partial 


differential equation, 


2 2 -- E 


Ba 21M^ Ky, - 0. απο 


This equation is nonlinear, of mixed elliptic-hyperoolic 
type, and contains variable coefficients. However, because 
flutter analysis is, generally, concerned with the stability 
of small perturbations about a steady flow, it is appropriate 
to assume that the oscillatory nonsteady component is small 
compared to the steady flow potential. Thus, the governing 


equation becomes 


(AMÓ Yo + yy — MCYFLI LO yy + Fy Vy 


= 21M Ky, t MK^w - 0, πι. 


with the usual nonlinear steady flow equation also resulting. 
Equation (II-4), though linear, retains variable coefficients 
dependent upon the steady flow solution, the nature of which 
(subsonic or supersonic) determines its type (elliptic or 


hyperbolic). 


B. REDUCTION OF THE GOVERNING EQUATION 
First, the reference Mach number shall be unity, reducing 


rus lon (11-4) to 
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byy = NED + Buggy) - 24K09z + K%9 = 0, (11-5) 


Second, the approach taken by Teipel [Ref. 30] and 
Hosokawa [Ref. 11] is utilized. That is, the nonconstant 
coefficients on Yyy and y. are approximated as follows. 
Pirst, the steady flow velocity perturbation is assumed to 
s Small in the X direction, implying that the airfoil or 


blade is of small thickness, 
dy -. 9. (II-6) 


second, the steady flow fluid acceleration over the airfoil 
or blade is assumed to be a constant (parabolic are airfoils 
exhibit nearly this characteristic), which allows the intro- 


duction of the Oswatitsch parameter 
NE (y+1) byy "τα. 20% (11 - 1) 


The Oswatitsch parameter was originated by Oswatitsch and 
Keune [Ref. 20] for steady flow, incidentally. Thus, the 


governing equation has, finally, become 


(uy 7 O21). * K^y = 0. πι ο 


Note that this equation is of the parabolie type. 
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eee BOUNDARY CONDITIONS 
1. Flow Tangency 


The usual surface boundary condition is 


DF 


DT = ο. on HKI T), (II-9) 


where F(X,Y,T) describes the surface of a body moving as a 
function of time. The boundary condition requires flow tan- 
gency, that is the total derivative of F is zero. Then, con- 
sidering a thin airfoil restricted to vanishingly small ampli- 
ὃν of plane oscillations (recall the small nonsteady 
potential argument), its upper and lower surfaces can be 


described by 


F(X,Y,T) = = H (X,T), 
(11-10) 


PQ(X,Y,T) - Υ - H,OGT), 


where 


EEX T) = th, (X) + R. P. [h, G0 e 777] i (11-11) 


u 


With the definition of $(X,Y,T) from equation (II-2) 
and H(X,Y,T) from above, the linearized flow tangency 


condition becomes 
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py(X,0t) = 3h(X) , ΤΙ (MASON) 
2. Zero Upstream Disturbance 
The governing equation nas become parabolic as shown 


by equation (II-8), requiring, upstream of the leading edges, 


W(X,Y) = 0. (11-13) 
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III. FORMULATION OF SOLUTIONS 


A. SOLUTION THROUGH THE LAPLACE TRANSFORMATION 


l. Cascade Geometry and Boundary Conditions 


po _— 2 ———À 


U 
m x 
E... 
T 
Eu m A e. Do E 
---ᾱ--ἡ -----.-.5ο 1 


Fig. 1. Unstaggered Cascade. 


Consider a pair of the parabolic arc airfoils or 
blades in an infinite unstaggered cascade (row) immersed 
within a transonic flow field. The blade chord is 2b, inter- 
blade spacing is P, the lower blade is the reference blade 
(variables associated with it are subscripted, zero), and 
the upper is referred to as the adjacent blade (subscript, 
one). The coordinate reference is fixed to the reference 
blade midchord for consistency with Teipel [Ref. 30] with 
whose work comparisons were made. Pitch and plunge ocsilla- 
tions can be described, pitch axis at a, magnitude of pitch 


Qo > of plunge ho. The adjacent blade performs a harmonic 
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piteh/plunge oscillation arbitrarily phased with the oscilla- 
tion of the reference blade. 
Two n coordinates are specified, with "ng and n4 


related such that 


t P, IE) 


> 


> 51 
ΠῚ 
n 
LL m - 0,1, 
(ED 
E 
Ü “` = > 
- 
ome pee 
ES 
1 b 


Note that the linear nondimensionalizing parameter 
Will be the blade semichord for all variables in this analysis. 


Thus, for example, the reduced frequency 


κ- bb 6 (Ma) 


el 





From considerations in Section II, the boundary con- 
ditions can be expressed, for the reference and adjacent 


blades, as 


3h 
= RE u 2 
Yy (X,Y 50) = vo (X) = . - ikh,, : ome 
and 
Y, (X,Y =p) 72 v, (GO = ev (x) Con } 
MO 1 0 > 


where G is the interblade phase angle. 
Then, for pitch and plunge of small magnitude 
(III-6) 


ΣΣ == agL 1*iK(X-X ie Ka 


0 m 


For computations for the steady flow case, a 
parabolic are airfoil is assumed (Y = mem see Ref. 30) 
and 


vo X) = 2 2TX. μμ ο 


EM SoIwmpronveot the Oscillatory Sonic Flow Equation 


The governing equation from Section II is 


Vyy ~ Wy + Kev = 0, (111-8) 
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where 
OA (MARIE O) 


For utilization of the Laplace Transformation define 


anew streamwise coordinate 
Von δις], Cpu 1ο) 


Denoting the Laplace transformation of a function 


with respect to X4 with an overbar, i.e., 
L e -pX, 
f(p) = Lf (X,)] = fe £(X,) aX. , (GET) 
0 
transform equation (II-8) to 
2 
Wy ντο, (III-12) 
where 
u = Qp - KÊ, πο... 


The boundary conditions become 


Y, (p,0) = V,(p), ο 


> 





by (P50) = e!?v (p). 15, 


The solution to equation ΠΕ 1) ουο οσους to 


MTS and 15) is 


πος... 
ΠΕ 2 e sosh(ut) - esehlult-p)] (MATTON) 


The hyperbolic functions are expanded, as follows, 


to allow a term by term inversion 


ο” «Ώμρ πμ. 


e = 28 "P y , 


T m EE 
sinh(yup) (eHO _ QT HP) n=0 
then 


( e uL-Y* (2n*1)p] + e uL Y* C2n*1)pJ4, ae) 


cosh(uY) . = 
sinh (yup) = 
and 


A παρ (III-19) 


cosh[u(Y-0)] . = - 
o) O 
HEUS. 


= -ui-Y-t ar - + +o) 
(¥,(p)e* Le ul-Y+2notp) E neam 024 
0 
CUTIES 210) 


00 
2 


= 1 
y(p,Y) = = 
Hn 


C7 (p) [e7 HC- Y+2np+2p) , ¿-u(*Y+2np)j), 
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For inversion, recall the formula 


co “PX, -o/lx, 


f e 
en Be QX. = Jr τ .π amnem 
° Xx; i WP 





Utilizing the above formula and the convolution 


pheorem, recognizing that 


NES Όχι «8, ( P REO ) 


and reintroducing X after inversion, yields 


2 


X Sy (X-S) " n e a 
Mey) A , (ser e: I 
ynQ -1 X4X-S8  n-0 
| Q(2np* o*tY)?  Q(2np*2p*Y)? _Q (2np- y) 
H (x3 — 1 (X-S) 
A EE v, (5) Le Kee "n x= 5) Jas. 
πιο) 


Then, for the potential on the reference blade 


Upper surface, 


k^ m Q(2np*o)* 
ποσα EN han e 9 τ » {v,(S)e*° [2e ΕΘ ] 
ViQ -1 O 
_Q(2npt20)* -QGng)- 
- v,(s) Le j 0.4.7 G-S qas. (III-2H) 





This equation was programmed for digital calculation (see 
Appendix A and Section IV-E). 
3. Modification for Sonic Free Jet Problem 

The sonic free jet is of interest as an extreme case 
pend tunnel wall porosity. That is, the porosity is so 
great that the wall has vanished and the airfoil is in a 
sonic jet constrained at its boundaries not by a zero 
velocity condition but zero pressure differential. 

The boundary condition on the reference blade remains 
the same as before, equation (III-4). The adjacent blade 15 
replaced by the edge of the sonic jet upon which the following 


pressure perturbation condition is imposed 
3Ky(X,Y=p) + y, (X,Y=p) = 0. 55; 


transforming and substituting as before, the second 


boundary condition becomes 
(iK + p)y(p,o) Ξ 0. (ión) 


ὑπ »ροΐη conditions, the solution to the 


transformed equation is 


2 


vo (Pp) „-2up.u! _ Zur, 
u o + 1 


v(p,Y) 





v (p) 
O curata on) 
nn) (111-21) 
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with 





1 s EU E IUE 
= = 26 Σ (ob) e 1-25) 
cosh(up) ES ch 5 
E Vg (P) E n, -u(-Ytanotp) -u[Y*2np] 
v(p,Y) = n 2 Kalle EMO -- A 209 
n=0 


The subsequent inversion is similar to that in the 
previous section. 

The result was programmed for digital computation and 
used as a check case against which the Collocation techniaue 
porous wall wind tunnel solution was compared. The Fourier 
Transformation method (see next section) was also applied to 
the free Jet problem and, by following a derivation similar 
to that above, it was found that an interblade phase angle 
of zero in the cascade problem correctly modeied the sonic 


free jet problem. 


B. SOLUTION THROUGH THE FOURIER TRANSFORMATION 


1. Cascade Geometry and Boundary Conditions 
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a. ee ee A DA E Es cs. eee a a 


Fig. 2. Unstaggered Cascade. 


The arrangement of the cascade is the same as the 
previous section, however, the coordinate system is changed 
to simplify the πο ὃς The origin of the sss 
system is now fixed to the reference blade leading edge. 
ade chord is c. 

The following nondimensional variables are defined 


for this phase of the study: 


52 
ΠΕ 
ye, 
(III-30) 
P 
zer 
_ 2 
_ 
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h M 
2 ο 
All variables will use c (vice b) for nondimension- 


allzation as appropriate, i.e., 
K. = — . (III-31) 
ine boundary conditions are the same as those expressed 


in the previous section, though their representation is 


somewhat different. 


Wy (X,Y=0) = w(X) , 

by (X,Y=h) AO ON (cOn 
and 

w(X) 9 [-1-1K (X-A,) Joy * iKjh, . (ue) 


ο ranstormation of the Sonic Flow Equation 
ειστε, Degine che equation for oscillatory 
sonic flow from Section TI and transform πη the 
following Fourier Transform, 


+o Min 
A Ae ax , ο. 


a 





Obtain 


Yyy + μ΄ψ ΟΠ; 0111-55) 


"a = — Ἱλν + «Κιν T Eu (III-36) 


This equation is satisfied by 


UY + P e iut 


v(v,Y) = Eje 1 


(III-37) 
Using equation (III-32) 
E 4σ-- 
Py (v,h) = e “py (y ,0) > 
Ehen : ë O) - p. [e 7 Hh E E 
or 
eic, iuh 
= 22 
i - i (111-38) 


Then, with 


Yy(v,0) > w(v) 5 
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E. and F. become 


1 T 

π. - πίν) - - οσα” μη. 

1 e g ο μη e iuh 

and 

een e la 

F; u al J 

iue Q iuh e, iuh 

«bmc 


mo). O) "(ice i Se BP) et HT MA e m 
3 a πάρκο EMG. τα ντ ο ο AR 
X: e iuh Άν ΤΊ 


] 
SO 


Because the pressures on the reference blade are of 


interest, set Y = 0. With 





ia, -iqa πο .. Τα 
cos a = 28 ee OS 
2 > 5 
and 
Y - CERE | 2 tan 
ορ = LE, Sin 2a = — — ; 
Jr tan Q L Ὁ ttan g 
then 
= E - uh 
T 2 W(v)r2(1-e Aes t 2(1lte lo tan! 2 
Ψ(ν,0) s 


ue 4tan 


St 





or 


)cot 2 " 








w(v) e 


(eias 1 
u 2 2 


)tanE2J, (IIT-40) 

It can be shown that equations (III-16) and (III-!0) 
Exx-guivalent and in this sense the two solutions are not 
independent. However, the inversions are substantially 
different. The following stems from the work of Hamamoto 
[Ref. 10] and, apparently, arises from observations of abstract 
eomplex variable theory, it differs from standard Fourier 
Transform inversion techniques and ieads to a very compu- 
tationally convenient result. A further method of inversion, 
utilizing the Theorem of Residues, was considered and is 
discussed in Section V-C. 

Eee Inversion 
In order to simplify the inversion of equation (III-40), 


Sy and Ty will be introduced and evaluated in a subsequent 


section. 
Define So and To as follows, 
1 μη uh ivX 
Sg X) = Dom E man Eo" e are (ITI-4]) 
and 
1 e i Hide a x 
Tg (X) Er y „cot Le aX. CO pana 2) 
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Also 


$ 
Sy+1 60) = E Sy n? dn, (III-43) 
and 
= ; 4y 
Tel) = CT) dn. (111-14) 
Using the convolution theorem and pesoonmzogne that 
w(X) = 0, y(X,Y)-20, X«0, (III-45) 
yields 


ig 49) 
i Έα) 50 (Φ) 1 ἀφ. (111-216) 


X 
TA O T4) * ( 
0 








me eerate Dy parts to obtain 


X 
ΤΟ; - -{w(O)[D-T, (X) i B-S. (X) ] iR ὶ Ww(X-0) [D T4 (0) 
0 
+ B-S,(6)] do} , (III-47) 
where 
io ic 
= - 1 _ wu T 
ra ==, p= —=—-. (III-!8) 
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4^. Calculation of Lift and Moment Coefficients 
This approximate solution of the sonic unsteady 
cascade probiem (III-47) was found to be quite efficient in 
its usage of computer time. Thus, it was extended, as follows, 
bo calculate pressure, lift, and moment coefficients in 


Euddtlion to velocity potential. 


Since 
and Ë 
1 | 
Cr = -2 "1 C5 (X,0) ο. (III-50) 
then 
- : 
Cr = l| E [y (X,0) + ik, w(X,0) J ax 
AL | 
= Uf f Vy (X, 0) ax + iK. Tero Xe OMEN (III-51) 
0 0 


With the previously defined y and w(X), and the 
blowing definition 
T 
E ο ο. ο ας O T GO (0 


“1K, [D-T,(1)+B-S.(1)]} , (III-52) 
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the lift is 


Cr = 4[w(1,0) + 1K,PF5(1)]. 


For the moment coefficient 


2 : 
2 E (XA, ) Cp(X,0) ax 


©2 
1 


i] η 
ο S Xp(X,0)dX + S Xy, (X,0)8X] + A.C... (III-54) 
0 0 Ze 1 


H 


Integrating by parts and eontinuing the notation 


started above 


Cy = =A-1K, 1F,(1)-P3(1)) + [9(2,0)-F,(1)7) + A, -C,. (LIL-55) 


2. Solutions for 5, 8nd Tw 
Because Sg (X) and Tg OO are difficult to evaluate 
directly, I(X) and J(X) will be defined and calculated. It 
will be seen that I(X) and J(X) are S. (X) and T. (X), respec- 


tively, from which all Sy and Tu can be obtained. 


Define 


to  ivX 
I(X) = NM J - tan = dy: CDS) 





tuu cem cp cv 
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TX 4X 


tan 22 = p — H 
2 _ m=0 (2m+1)°-X° ° 











with 
X h h 
SAA men, 
E tan Th I us Σ I 
2 ο... E 
2 2 ne 
tenen 
t9? iyX © 
IX) -54p { I; Ξ =] av. 
o ο ο. 2n-1,2 m^ 
(5) - ( o —- 
2 2 2 
h 
Substitute for p, using equation (III-35) in the 
bracket 
2 : 2 "gk E x | 
= = + + E . 
n=l -idAv eR_v K. d ¿¿n-1, 2 ne Tn (1-33 )(v-v_) 
Y 2 ne μη 
2 2 
E ει νζε-ιλι, (111-51) 
n h l T 
SO 
1 1 οο + co ¿Vx 9 
E E ’ τπππην' pene 





Now define 








5 to 1νΧ 
Ma ph | 
J (X) 3T m Sy cot “> dv. (III-59) 


With [Ref. 12] 


eo 





IO = I + 2χ΄ Σ —— 
τ 2 2 
m=l X - v 
obtain 
πο iux co 
j| 2 1 
J(X) = —- jg [+ E =] dv 
2m1 Be vh 1 n=2 ,u.2 (n-1)° - 
(5) - 5 
h 
With 
2 2 
πι 7 2 . 
De L E n, (III-60) 
then 
co + co LX 
»i D Al E u 


The integral occurring in equations (III-58 and 61) 


is, from Campbeli and Foster [Ref. 4], 


πο. lux İv X 


e Em £ 
Ao vO-V 7 dv = v [e - i (111-62) 
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Thence 


E İv X 
2 (1 — e ) 
T(X) = I c 
h(eK,-iA | m 
and 
oo A 
er Ze: x -ny Lj = — 
1 n=1 BD 
Now, since 
ϕ 
S (9) = J S. (n) dn , 
0 
then 
] E JL uh iv n 
ae - S, : 
S40) 2m i L y tan Ze Maul dn 
+ 00 τη 
ES eh n 
> n tan > y E dn! dv: 
Therefore, 
S,(9) = 1(9) 
Similarly 
T, (9) = J(9). 


Hl 


πιο) 


(III-64) 


(III-65) 


πα ος) 


(III-67) 





μας, 


2 co 1 iv X 
S. (X) = Seem 2 m == MS ] 9 Μπιτ ος) 
1 Dorn I 
and 
X 
S,(X) = f S, (n) an, (III-69) 
0 
5ο 
co x iv 
2 L n 
S AX) = — ; — / [l - e ilm 
e h eK, Se D 
co iv X 
= 2 5 τ ES a n 
Rao + X o* ge ELS 
i n=] n n 


(III-70) 


ο ο. and likewise for Tr 


C. SOLUTION THROUGH COLLOCATION TECHNIQUE 

Utilizing the same geometry, nondimensionalization, and 
Dommgsey conditions as for une first portion of Section III, 
the boundary value problem is solved by assuming the potential 


function is a summation, 


W(X,Y) = $6 (X, Y) * $505 Y) + $2 Y) * 0,1), G 
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where $6 and N are the isolated blade potential solutions 
and bo and $1» are the interference potentials for the 
reference and adjacent blades. 
Such an approach was proposed by Gorelov [Ref. 8] for 
sub/supersonic flows. See also Samoylovich [Ref. 21]. 
Utilizing the boundary conditions and the known source 


distributions 

0 _ E 

$5 y (6,0) n vq UO > m = 031 5 (III-72) 
obtain two equations 


=-p) = 0 


ISS) 
Reference Blade, 


0 
Y = Y =— + Y 
ος om ο a a 


0 = ? ( = = = 
Poy (EY en SUS) Qy (X.Y 0) = 0 


Adjacent Blade. (III-7) 


The isolated blade solution is already known from Teipel 
Mer, 301 so 
Ke 1 
X | Expl=7(X-8) - 7x25)] ° 


0 Ν 
φρίΧ»Υρ) EN / V S) 


— . f vg S, (111-15) 
Ὕπῶ -1 ΝΧ-5' 


46 





and 
2 
2 QY 
1 X Exp[ A (X-5) - i ] 


$2 (X, Y.) > (s) X-3)- as. EMI 
πα -- ARES | 


It is expected that interference potentials will have 


the same form but contain unknown source distributions u 


0 
and u, Thus, 
2 2 
᾽ X Exp[(X-S) - N cm ] as 
$g (X, Yg) Sa if Uys) SS >. 
MA => SS (III-77) 
and 
2 2 
K Q 1 
X Exp[—-(X-S) - σσ dS 
$ (X,Y) = o if u (s) _ ES (nuc 
-hQ -1 AX-S' 


HN τηε fundamental point of this solution to approxi- 
mate the unknown and known source distributions by series, 
the coefficients of which are to be determined by satisfying 
emat ions (111-73 and 74) at a sufficient number of points 
eche blades. 


Elder [Ref. 7] utilized a power series, 


Bor - Pow xt, (III-T9) 


ai 





and 


eoo 
vb) = E Ὁ X Ππ 0.1. (III-80) 


n=0 
Elder obtained the coefficients of these series. Expanding 
upon his methods, this author obtained potentials. For small 
Boubers of collocation points (L«10), 170, K«1, agreement 
with the potential calculated by the other methods was 
generally good. Numerical difficulties above those parameter 
Panges and, especially, for A=0 led to the introduction of 


the 'step-wise-linear' series approximation 


L-1 
= XY 4 ` = n - y 
ug X σῷ ο Y + I Vg X X, )H(X X, ) : (III-81) 
and 
L-1 
== ° + = = = 
ELO) 90 + ol X x OS X )H(X X) 3 ποτ ον 
where 
Ü 2 «0 
H(z) = (IFI 655) 
IA ο 


The first two terms were necessary to allow for initial 
value and slope of the source distributions (subsequent 
coefficients represent changes in that slope at the selected 


points, XN). Further, they allow simple representations for 
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the isolated blade source distribution coefficients. For a 


eech oscillation about the midchorá, 


where 


and 


where 


vg OX - 900 + 8 


X 


one ws 
ρε. (III-84) 
E 

9,, 7 - iKe^?, (ES) 


Inserting these step-wise-linear approximations into 


Exusvrons (LII-75, 76, 7T, and 79) and substituting the 


resulting expressions into the boundary condition equations 


for each blade (equations III-73 and 74) obtain 


E 





X E LE x s)-8 nn 
JS v f eg AS Gone 
ΗΠ, ο ——- 


ρ 
w 7 


2 2 
K Q P 
P8 X Εχρὶ-σίΧ-8)}-ῃ Τπτοτ. = 
r 2 T a J S Q OSS d ME UID ο (X=) 

zu (x-5) 3/2 > On n n 


X (S-X )H(S-X ) n 
m = ES expt (-8)- πο | as} 


ο [6 
Eo v 


x xp R(x. Ses ] 
= 0 [a Y t 


X _ Exp [É O- . E 

IO p ds 
s —T - 

DONE (x-8) V? J 


Reference Blade, (III-86) 


[Κ΄ -8 οὔ ] 

X Exp| —I(X-S)-r ais 

g ye Q Ἡ (X-S) UA 
En 2 yn "oo .^ t X 





— 7 απ. 
G. 11 
X exp Es) ] as co ' 
4 £8 v f S ὧν ες = E Lv, (X-X_)H(X-X_) 
27407 Ol -1 (X- s)3 2 deer Tn n n 


+ 


sm: X (S-X,)H(S-X ) 
A al ΚΞ exp O- 3-4 Aalen) 


rs Q ρ΄ 


2 IT o0 3 — — (X-s) Y? 
2 2 
á K Q p 
X Exp[—-(X-S)- CO ] as 
01 nd S E z — 
u 7 METEO 
Adjacent Blade. 


716 
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These equations are satisfied at each coilocation point 
to obtain the unknown coefficients. The coefficients can 


then be used in the following equation for the potential, 


1 X exp (x. S)] 
= ea {L859 f Yoo / Ze es 
AnQ - (X-5) 


¥(X,¥)=0) 
0 T 


Ke 
X ExpUg (X75) ] 
1 


172 4» 


αυ + yy. ] wenn 
01 01 -1 (X-S) 


X (S- X SESS Aa ) 
HO Exp [É (x- Sa as } 
(X-S) 


p? 
X Exp ex e 2. 
f T DES) 


IO 10 (X= 5) 172 m 


ὝπΩ τη 


X Exp E S)- - Posy] , 


+ [0..+v..]1 / 8 E Ça 
MED ᾿ uer 


co X (S- X, )H(S- x, ) Ke Q p^ 
τω. v jf expli (X-S)- deus 
7 SE Y 
(III-89) 
IE NAL UATION OF THE INTEGRALS OCCURRING IN EQUATIONS 
(III-24,86,87 AND 88) 
The integrals encountered are of the form 
X N Ke ap“ 3 
= r πως. = ES ΤΚ. ° 
In E c SS Expl- (X =) 0-57} do, (III-09) 


E 





nr oqnin the substitution, U = 


X-5, separating I 


N 


into its real and imaginary parts, and simplifying as did 


Elder, obtain 


E _ p (x-u)N 
NR 0 ym/2 
and 
X+1 (κυ) 


ΝΤ d m/2 


where 


1° +4KS 


Nenas necessary to evaluate six integrals 


X+1 
me - y 
N° ο 


Mene j = 1,2,3. 


ΕΧΡ[λΑ0-ΠΊ cos[2KAU+Ë] dU, 


Exp[ AAU-7] sin[2KAUt=] dU, 


CL αν οπ. 
ES _ Kp? 
5 E m 1:3 I Ξ T 


COS 


E 
sinl2KAu+T] au, 


Exp[ AAU-Z] 


(III-90) 


(III-91) 


(III-92) 


(111-93) 


(III-94) 


These are four cases of treatment of these integrals: 


l. The Steady Case 


Setting K = 


O reduces the integrals to a relatively 


simple case which may be evaluated through the use of the 


error function or by numerous numerical techniques, one of 
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whieh is utilized in the routines described in the 
Appendices. 
2. The Oscillatory Case with X = 0 
In this case the integrals are of the form, 


χει , 
"ιτ uJ 2/2 COS [2KAU+Ž] dU, επ») 


e 0 sin 
requiring a distasteful integration near the lower limit for 
J « 5/2, due to the rapid oscillations and increasing magni- 
tude of the integrand there. 
homevOld=eGdihittienr oles ab the Lower limit, the range 


pi ntegratíon is dissected, 1.€., 


Dok XA 


J=5/2 908 yE s j-5/2 cos E 
M U 33 pp C2KAU a au s U f ae au 
x JU 
j-5/2 cos E EC 
T n U 33 n L2KAU π. du, (III-96) 


XA 


The limit XA will be chosen so as to allow the second 
integral to be evaluated numerically while not greatly af- 
meeting accuracy. The first integral, hereafter the 'correc- 
noxuntegpal'. wilXi be approximated in the manner that 
follows. 


Recall 


cos [2KAU+=] = cos[2KAU] cos [=] = ela eA sin[=] ; 
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and 


sin[2KAU+FJ = sin[2KAU] cos [zJ + cos[2KAU] sin[Z] 


and, the truncated series expansion 


2 
eos[2KAU] 2 1 - E E 
and 
| Ν (2KAU) > 
sin[2KAU] 2 2KAU - LAI 


Then, for 2KAU sufficiently small 


ΧΑ ΧΑ 


J-5/2 cos E - j-5/2 cosrfE 
ad U sin 2RAU+TIAU N U a au 
DER 9 e. 3 XA + 
£ j-3/2 sinrE (2KA) j-1/2 cos;-E 
+ (2KA) E U cost au un τ U mel 
XA . š 
(2KA) j*1/2 sin,E 
Since integrals of the form 
XA 
m/2 COSE -98 
fou sig au (III-98) 


0 
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may be expressed in terms of the Fresnel Integrals, the inte- 


grals in (III-97) can be calculated after integrating by 


parts 
y4-5/2 cos EE? j-3/2 cos. E 
P" Ping ldU 2j- agg sin XA- 
A 3/2 
2 U MESE 
+ 2E "d ου Με E aus (111-99) 


where, j > 3/2. From Elder [Ref. 7] for j = 1 


XA 


f v? cosc Jau - Ns EU)» ευη, 
0 
where 
ES cos 
E XA MI» dt 
NM LS ο. Πτα) 
0 Se 
Then, for j = 2 

XA 3/2 sin-E- 

= U; 
IU 1/2 cos [Sau = 2-XA E cost + 2°E p “u” cos aA ie 
0 0 


(111-102) 


Thus, the correction integrals can be evaluated sequentially 


on increasing j. 
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ΠΕ General Case, K770 and AFO 


ln this case, as the lower limit is approached, the 
exponential function in the integrand is depended upon to 
ave the integrand magnitude to zero more rapidly than it 
is magnified by the negative power of U. ποια Na ORG Lon 
Of the integral near zero ls neglected by specifying a non- 
zero lower limit for which a negligibly small error 45 pro- 
emecd while still allowing numerical integration to be 
terized, 

H. The Case for P=0, j=2, K¥0 and 470 

ο ας Case, another approach for that portion of 
the integral near zero must be taken. Due to the absence 
of the -B/U in the exponential function, that portion cannot 
be neglected as in part 3., above. 

The range of integration is again dissected, a 
‘correction integral' is defined, and XA is selected simi- 
fey tO part 2., above. The ‘correction integral’ integrand 
is expanded using the series expansions for tne functions 
within and higher order terms in XA neglected. The resulting 
series is integrated term by term and programmed for correcting 
the numerically calculated portion of the integral (see 
Subroutine ZSUM in the Appendices). 

E. APPLICATION AND COMPARISON OF SOLUTIONS TO SEVERAL 
CASCADE PROBLEMS 
In order to gain some confidence in the Collocation method 


solution before extending it to the porous walled sonic wind 
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tunnel problem, the solutions obtained for sonic cascade 
interference in the previous sections were programmed for 
digital computation (see Appendices A-C). They were applied 
to cascade problems differing in the parameters of interblade 
phase angle, reduced frequency, Oswatitsch parameter, and 
interblade distance with pitch oscillations about the cascade 
pesseomidchord position. Many of those results are shown in 
this section to illustrate the ranges of parameter values for 
which agreement between solutions can or cannot be expected. 
Because it is the object of this section to illustrate the 
basic traits of the solutions, especially those of the Collo- 
cation method solution, agreement will be considered to be 

a 5% or so difference. 

The complex velocity potentials calculated by the three 
methods are presented in three sections. First, reduced 
frequency, Oswatitsch parameter and interblade distance were 
fixed, the interblade phase angle varied. Second, interblade 
distance and phase angle were fixed and the remaining parameters 
varied to illustrate their interrelationship to the agreement 
between the Collocation and other methods. Third, the inter- 
blade phase angle was fixed, the Oswatitsch parameter set to 
zero (the governing equation is thereby linearized), and the 
reduced frequency and interblade distance varied to illustrate 
Moiewecwificulties the Collocation method had satisfying the 


linearized transonic potential equation. 
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a no I] S shown indicate that the Collocation method, 
as applied here, agrees reasonably well for lower frequencies 
and larger Oswatitseh parameter with the other methods. For 
larger frequency and/or small Oswatitsch parameter, the step- 
wise-linear series approximations to the then oscillatory 
source distributions become oscillations whose amplitude may 
grow substantially from the blade leading edge aft. For 
cases in which the amplitude growth is large, agreement be- 
tween the potential calculations of the Collocation method and 
the others deteriorates. Furthermore, increasing the number 
of Collocation points fails to improve agreement. Finally, 
it should be noted that when agreement was poor, the source 
distribution series approximation yielded an oscillatory 
result whose amplitude growth was large. 

It has been observed that a step-wise-linear series 
approximation can exhibit this characteristic when numerical 
errors of sufficient magnitude are present. Further, the 
effect of increasing the number of Collocation points (and 
coefficients) is to increase the growth once it was present, 
as observed above, Numerical errors of substantial magnitude 
can be traced to the calculation of the integrals (their 
values are highly damped oscillations along the chord) in 
the Collocation method solution above. 

Maire bo produce accurate potentials when an oscil- 
latory series approximation was obtained can be attributed 


to the large magnitude of the coefficients near the trailing 
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edge. The resulting potential calculations will then involve 
the difference of large numbers. 
ie tnberbtade Phase Angle Variation 

Table ieamews real and imaginary potentials for 
Ieee. tO, A = .10, p = 1.0 and o = T, ı/2, and O radians 
(recall, o = m corresponds to an out-of-phase oscillation of 
the blades) for the Collocation, Fourier Transform, and La- 
ο. ο Transform Solutions. Further, the Collocation and 
Kere methods are shown for K = .10, A = .50, 0 = 1.0, 
andio = 0, 

Agreement is seen te be fairly good for dg > 0. At 
See, Coe Coliocation method fails to agree with the others 
οἱ about X= -.10, unable to fit small amplitude imaginary 
part's oscillation. An increase in thickness (increased 
Oswatitsch parameter) improves the agreement slightly. 

2. Frequency and Thickness Variation 

Figure 3 shows the results of potential calculations 
versus position along the cascade blades described in Section 
III-A (-1 < X « 1), for Collocation (step-wise-linear series 
approximation), Fourier Transform, and two cases each of the 
Laplace Transform and Collocation (power series) methods. 
Interblade distance, rho, was fixed at one and interblade 
phase angle, o, at pi radians (out-of-phase oscillation 
produces greatest interference). Reduced frequency, K, was 
EE ua adues .01, 1$. .3, .5, 1.0, and 1.5. The Oswatitsch 


Eu vcr Look velues..05, .1, .2, and .5. For all but 
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TABLE I 


COMPARISON OF POTENTIAL (REAL/IMAGINARY) CALCULATIONS 
OF COLLOCATION, FOURIER, AND LAPLACE SOLUTIONS FOR 
VARIOUS PHASE ANGLES 


X COLLOCATION FOURIER LAPLACE 


EEESO TO. 1 = 0.10, p = 1.00 


For o = 3.10 rad 


=.9 ο ο - 96 0.486/ -0.854 0.489/ -0.854 
-.7 1.048, -2.5ΤΗ 1.146/ -2.525 AO p OY 
5 AL. aL SIN Baie ea -u96 15» 
-.3 2.426/ -5.937 2.530/ -5.860 AO O 1 
-.1 a = 1.007 pP MA 3.264/ -7.508 
De BOO 29.266 O ΞΟΠ οἱ ο τς 
0% 3 Wowie = 1.0 - 920 4. 768/-10. 834 A 8 10 
9. 5 5.459/-12.570 5.558/-12.486 6506027212 ου 
ουσ 6.271/-14.220 ο πο τη το. 6.392/-14.110 
0.9 7. 104/-15.870 προ 15 27080 τοις 0 
For o = 1.60 rad 
- ο ο =0. 448 ο οσο 6. 16 0.094/ -0.404 
m7 -0.497/ -1.574 -0,425/ -1.605 -0.423/ -1.604 
5 -1.005/ -2.783 5 Το -0.918/ -2.776 
=.3 -1.503/ -3.927 = πο 23.9816 -1.403/ -3.956 
- 1 -1.960/ -5.156 -1.883/ -5.144 -1.872/ -5.142 
i. ος -ὂ οὓς ΠΠ 5» E SS 
0.3 -2.865/ -7.589 ΠΗ ΤΙ, τ; 
0.5 222107 28.5870 οι. ου ΗΠ 2392715 
047 ECL ale mL 556 za AOT 
0.9 -4,187/-11.000 nel] 122.206 -1,039/-11.210 
For o = 0.0 rad 
-.9 υπ υ 5 =o. 009 0.545/ -0.009 
-.7 0202 720.029 0.498/ -0.03% 0.500“ 0 055 
-,5 ο που  -ο ϱ-τ 0.499/ -0.025 501,7 20,025 
-.3 IO 01 002 0.499/ -0.015 50/20 016 
— 0.490/ -0.033 ORO, FEO ος 0.501/ -0.006 
01 "υπ ο ο 0.005 0250175592005 
0.3 0.466/ -0.052 029077092015 0.501/ 0.014 
0.5 0.563/ 0.129 0.499/ 0.024 0.501/7 0.024 
Qu T 0.393/ -0.128 Doo ο ου ο ο. ο 
0.9 OESO 0.296 0.499/ 0.044 ooo 5 ος 
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K = 


0. 


For G 


O@o @o I1 1 


WU Uy PP NU - NO 


TABLE 


COPROSATIÓN 
LOA T ο δρ Ξ 1.00 
= 0.0 rad 
0.434/ -0.088 
0.509/ -0.042 
0.489/ -0.024 
0.515Z -0.023 
0.485/ 0.001 
0.520/ -0.009 
0.480/ 0.030 
ο ου 0.000 
ο "το 0.063 
ο ο. 5 ος 
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500/ -0.039 
499/ -0.027 
MA OO 
HS OOO 
499/ 0.003 
499/ 0.013 
110072202025 
199/ 0.033 
499/ 0.043 
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Fig. 3-c, in which twenty were used, ten potentials were 
calculated along the blade and those positions were used as 
node points for the Collocation routines. 

Fourier results were machine plotted as a solid curve 
connecting individual points, which can result in a jagged 
curve for some parameter combinations because of plotter 
imitations. Collocation (step-wise-linear) results are 
plotted as crosses. Laplace Transform calculations were 
made for K = .01, À = .05 (computer time usage = 6 minutes) 
and K = 1.0, A = .1 (usage = 40 minutes) and plotted as 
squares. Power series Collocation results were also plotted 
as squares for K = à = .1 (usage = 14 minutes) and K = .5, 

À = .1 (usage = 50 minutes). Typical times for similar cal- 
culations by Fourier routine are 2-5 seconds and step-wise- 
Iumear Collocation, 10-20 seconds. Interestingly, though not 
surprising considering the method by which integrals were 
pEueulsted, the Collocation results for 20 points were obtained 
in slightly less time than 10 points. 

Figure 3 is arranged in order of ascending frequency 
D Nw varying at each to show the characteristics of the Collo- 
eacıon method solution, Generally good agreement in the cai- 
culation of the interference potentials for all four methods 
used is shown. Because Fourier and Laplace agree very well 
at the points shown and for many others not shown, the Fourier 
solution was concluded to be a good standard. It was further 


coneluded that increasing frequency adversely affects the 
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meeuracy of Collocation results and increasing thickness, as 
represented by tne Oswatitsch parameter, improves it. Further- 
more, increasing Collocation points does not necessarily 
improve accuracy as illustrated by Fig. 3-c compared to 
Ea. 3=D. 

It was observed that in those cases for which rela- 
tively poor potential agreement between the Collocation and 
the other methods was obtained, the step-wise-linear series 
source distribution approximation was an oscillation whose 
amplitude grew toward the trailing edge. The coefficients 
EE nose series, because they represent the changes in the 
ο ο Of the approximation, can illustrate the point. For 
example, for K = .3, A = .05, which shows good agreement, the 
first and last coefficients are (real/imaginary) -.7280/1.760 
and 1.488/-.0269, respectively. For K = .5, » = .05, poorer 
agreement, with -.1357/1.879 and -1.664/-2.942. For K = 1.0, 
A = .1, very poor agreement, with 3.549/1.402 and -164.4/-392.7. 
EE sand À = .1 and 20 collocation points, poor agree- 
ment, with -.4969/3.505 and 1714.0/-2466.0. The same result 
is seen in the reduction in accuracy shown for o = 0.0 in 


the previous section. 
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Yscale = 2 units/inch Yscale = 10 units/inch 


Real Part Imaginary Part 
Fig. 3-a. Complex Velocity Potential Comparisons for o=T 


and p=1.0 with Varying Frequency and Oswatitsch Parameter. 
L-Noedty Potential vs. X (1 unit/inch). 
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Real Part Imaginary Part 
Fig. 3b. Complex Velocity Potential Comparisons for 0=T 


and p21.0 with Varying Frequency and Oswatitseh Parameter. 
πι; Potential vs. X (1 unit/inch). 
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luc Complex Velocity Potential Comparisons for o-T 


and p=1.0 with Varying Frequency and Oswatitsch Parameter. 
Eco "Potential vs. X (1 unit/inch). 
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Wee as ποπ 


Real Part Imaginary Part 
' Fig. 3d. Complex Velocity Potential Comparisons for στΞτ 


and p=1.0 with Varying Frequency and Oswatitsch Parameter. 
ποῦν Potential vs. X (1 wnit/inch). 
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Real Part Imaginary Part 


Fig. 3e. Complex Velocity Potential Comparisons for G=T 
and p=1.0 with Varying Frequency and Oswatitsch Parameter. 
Welocity Potential vs. X (1 unit/inch). 
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Bes 3. Complex Velocity Potential Comparisons for o=" 
and p=1.0 with Varying Frequency and Oswatitsch Parameter. 
emocie Potential vs. X (1 unit/inch). 
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3. Oswatitsch Parameter = 0.0 

Fig. 4 is a compilation of selected interference 
Peg@ential calculations for A = 0.0. Because the accuracy of 
the correction integral technique used in the Collocation 
method (and Laplace Transform method) (see Appendix A-II), 
memmects accuracy to frequency and interblade spacing, various 
interblade distances were chosen to illustrate the effect of 
frequency on the Collocation method solution. 

Again, the values of potential are plotted against 
100, X, along the cascade blade (-1 « X « 1), Fourier 
Transform results are plotted as a solid line connecting 
points, Collocation (step-wise-linear series) results are 
plotted as crosses (except at p = 12, where Teipel results 
for an isolated wing are so plotted), and Laplace Transform 
results are plotted as squares. The number of collocation 
points varied from ten to twenty, as can be seen from the 
Pots, 

Fig. 4-a shows results for interblade distance rho 
equal to one. For K = .01, the Collocation results were So 
bad the real parts would not plot on the graph. For K = .1, 
the real parts of the potentials agreed fairly well. The 
bottom plots show the effect of increasing the number of 
points. It is seen again that the Laplace and Fourier 
Transform solutions agree rather well, as indicated in the 


few results shown and others not shown. 
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Fig. 4-b and c show results for p = H (p is increased 
EN nat frequency can be increased while still attempting to 
maintain accuracy in the Collocation routine, see Section 
III-D and Appendix A-II). Agreement is only fair and the 
conelusion must be drawn that the derivatives of potentials 
needed for calculations of pressure coefficients, lifts and 
moments are likely to be very questionable. As in the pre- 
tomes part of this section, large growth of the source dis- 
tribution approximation amplitude would signal poorer 
agreement for the Collocation method. 

The bottom graphs of Fig. Ta are for p = 12.0. The 
crosses depíct potentials calculated by a routine following 
the work of Teipel for an isolated wing. The convergence of 
Fourier Transform solution results into isolated wing values 
for increasing interblade distance will be discussed later. 
Here, these data further serve to illustrate the applicability 


Of the fourier Transformation solution. 


T1 





pn 






Am 


"σαι = 1 unit/inch Yseale = 100 units/inch 





Meeale = 1 unit/inch Yscale = l0mm»its/ineh 





Yscale = 1 unit/inch μοι To uJ EE US en 
Real Pant Imasınary ame 
Fig. 4a. Complex Velocity Potential Comparisons for 0=T 


and A= 0. Varying Frequency and Interbiace Distance. 
εἴν Potential vs. X (1 unit/inch). 
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Fig. 4b. Complex Velocity Potential Comparisons for g=T 


and A=0, Varying Frequency and Interblade Distance. 
ce lty Potential vs. X (1 unit/inch). 
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Fig. Ue. Complex Velocity Potential Comparisons for o=T 


and A=0, Varying Frequency and Interblade Distance. 
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Fig. ad. Complex Velocity Potential Comparisons for σΞπ 
and A=0, Varying Frequency and Interblade Distance. 
Pevocity Potential vs X (1 unit/inch). 
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IV.  TRANSONIC CASCADE AND BIPLANE STABILITY 


A. GENERAL 

in the previous sections, the Fourier Transform solution 
has emerged as a standard. The method has been made flexible 
by the inclusion of the Oswatitsch parameter and variable 
interblade phase angle. It has been compared to Laplace 
ieemostorm, Collocation method, and Teipel [Ref. 30] solutions 
for many different sets of parameters and has produced good 
RESULTS. 

The Fourier Transform solution has been applied to the 
problem of transonic, unstaggered cascade aerodynamic stabili- 
ty. Selected results are included to show the effects of fre- 
τιν, Dlade thickness, interblade distance, and interblade 
phase angle. 

The biplane problem was solved by utilizing the Fourier 
Sewe Oon tor the interference caused moment upcn the upper 
surface of the lower blade and adding the Teipel [Ref. 30] 
solution for the isolated wing moment on the lower surface 
sr Chat blade. 

Isolated wing lift and moment results published by Nelson 
and Berman [Ref. 17], Teipel [Ref. 30], and Hamamoto [Ref. 10] 
were utilized for comparisons. Barton [Ref. 2] and Rott 
[Ref. 21] were also checked. The Fourier Transform solution 
reprcduced all of the aforementioned results with reasonable 
accuracy, even when the interblade distance was about 20 


(necessary for comparisons when A = 0.0). 


76 





DU LRANSONIC, UNSTAGGERED CASCADE STABILITY 
Calculations of stability boundaries were made for certain 
Euperblade phase angle, c, blade thickness/chord ratio, τ. 
and interblade distance, p. The moment coefficients were 
then calculated for varying frequency to find the frequency 
at which the moment's phase angle was “ 180 or 0, for each 


pivot location. Blade thickness was converted to Oswatitsch 


parameter as per Teipel [Ref. 30], 
A = 1.559 [(y+1)1]° 3 (1V-1) 


The boundaries are plotted so that stable ranges lie above 
the curves which represent the boundaries of stability for 
that p. 

It is seen that for p = 1, the boundary frequency lies 
υπ curves for forward pivot locations. Further, fre- 
quency, interblade distance, and blade thickness increases 
were stabilizing as were interblade phase angle reductions 
and rearward pivot shifts. It should be noted that the effect 
of interference was destabilizing as indicated by the strong 


effect of decreasing interblade distance, P. 
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fee TRANSONIC BIPLANE STABILITY 

A single result (Fig. 5-f) was included to show the effect 
of interference on the Biplane. Because the blades have iso- 
eea blade pressures on their 'outer' surfaces, the effect 
of interference is much less evident than observed in cascade 
flows. This can be seen by comparing the cascade and biplane 
curves for t = .01, and o = m(Figs. 5-c and 5-f). For close 
interblade spacing, where interference is most prevalent, the 
biplane stability curves drop to lower frequencies at more 


M ward pivot locations, indicating greater stability. 
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Peo eo LON OM THE COLLOCATION SOLUTION TO THE 
POROUS WALLED TRANSONIC WIND TUNNEL 


A. PROBLEM FORMULATION 
1. Geometry and Boundary Conditions 

The geometry considered here will differ in that the 
adjacent blade of the Collocation method solution for the 
sonic cascade is now a flat wall, the porosity of which may 
De varied. 

The boundary condition acting at the reference blade 
(or, now, the airfoil) is as before 


3h 
E = = u : 
0 (X,Y 70) = v. (X) = y + ikh, (pa pal) 


9 X 
From Drake [Ref. 6], the condition at the wall can 


be stated as 
De M o, pl ikw(X,¥=0) + y, (X,Y=p)] = 0. Cua) 


The condition of zero upstream disturbance holds as 
previously. Furthermore, non-dimensionalizations are the 
same. 

2. Solution of the Oscillating Sonic Flow Equation 
As in Section III, the source distributions shall be 


assumed to be series. However, because the 'adjacent blade' 
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is now a flat wall, its 'isolated blade' potential and 


ESPJIvatives vanish. Thus, 
EU ! 
U(X,Yg) = 69(X,YXg) + do X, Yo) + $ (X,Y). 52 
Now, from the boundary conditions and 


(SNC) cs sp Ὁἳσ Cee) 


0 
Poy m 


the airfoil and wall equations become 
$gy (X, Y970) + 9, y (X,Y,=-p) = 0, (V-4) 
'and 
φῦ Ot Yon o) + $gy CX, Yg7 0) + $4 y (X, Y. 70) + 
σηρί[Φρχίχ,Υρ-ρ) + $gy OG Ygno) * eq (%,Y,=0)] + 


1KLéo (X, Yo 7p) + 9,(%,Y,70) + 9,(X,Y,=0)3) = 0. 
(V-5) 


3. Derivation of Expressions for Derivatives of the 
Bowen & aa 


In Section 111-C the isolated and interference poten- 


tials are expressed by integrals of the appropriate source 
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Ame tributions (see equations IIl-71, 75, 76, 77, and 78). 


The wall boundary condition requires derivatives on the X 


variable of those potentials. Thus, the following derivatives 


must be evaluated, 


2E E [ET ex uS p“ J 
ao , Wo U T OS) 
9X XN 4x78 


PIN = dion (V-6) 


ορ n 


Me es, 2 P Ὁ 
I. "ο Y (X-8]- 48 


(V-7) 
9X XN X-S 


PIO = 


Eor n = 1, 


2 2 

K Q P 

X Expl=-(X-S)- - ] 

RES) as 
XN XS 


22 


I 3x 


2 
Q P 7 


2 
K 
X Exp[—(X-S) - + 
0 ο ο ο (1-8) 


2 ή ELA 


XN (X-S) 


Thus, derivatives are required of two integrals 


E 


DIM - 3x 


X 2 2 
K = q 
E (X-S) "Exp [^g (X-8) ~ t Ceres ls: Ic 5 


87 





Applying Leibniz! rule 


da 


q q 
— $ f(s,ajas = f = f(s,a)ds d Ced - f(p,a)SP. 
D D a da d 


DIM becomes 


X 
DIM- / -πίί.- 9) Ἑ PL i Se Jas 
XN : i des 


CE: 
~ spas) - $ nl 


To differentiate the integrand, let U = X-S, 


3g UU ExpL AU + B/UJ) = 


{my® 1 + K^ jm : Qp^ 
Q n 


2 
E 2 pl. 


There are two cases to consider, m = +1/2 that is 


- X-X 
2 n 2 n 4-5/2 
NEN -1/2 Qp U 
DI- * ο) U Paoi C NE 5 expLC] AU 
E 
X-X P 
n Exp CE A 
-ᾱ / U M gxp[c] au + Tt , — (V-10) 
0 


and 
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2 ni 2 n 
he -1/2 QP 23/2 
DIt Š T U ES EE ES P. U Ετος οὐ 
X-X 

n 2 0 

ο j υ-/πχργο]αυ ἐ Π/΄εχρ[- Es. (v-11) 

2 T 
where 
e Qr 

Q E ue 


The last term in equation (V-11) vanishes nicely for 
all P. As long as P > 0, the last term in (V-10) does as 


NEO: Pp = 0, however, care must be exercised. 








For P = 0, 
Ke me 1ο. Κ΄ 
DI-(Pa0) = =- | U Exp[—U JaU 
Q 0 Q 
χ-χ 
n 2 0 
MS -3/2 ArÁ l 
-- 5 Bd U Expl- UaU d -7;| 
U 
Since 
ax ax ax 
e ή A a e ET 
/ m Sia ed c m-l _m-l ü m-1 J n-] ns 
X X x 
then 
X-X ¡CAES 
= Expla U] N n 


n 2 
—3/2 K 


„2 
2 AA, Es U] 
m rU 
Y p „1/2 | 
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aus, 


x? 
Exp rn ] 


'' O) (V-12) 


B. RESULTS OF POROUS WALL TUNNEL CALCULATION 

The results of the preceding analytical work were pro- 
grammed (see Appendix-D) and pressure coefficients calculated. 
Comparisons were made with solid wall tunnel and sonic free 
Jet results obtained from the other methods. Experimental 
comparisons were made in the steady flow case using results 
Smeopreiter, Smith and Hyett [Ref. 27]. 

Figures 6,7, and 8 plot solid wall tunnel program results 
as crosses, sonic free jet 85 squares, and porous wall pro- 
"ΠΠ results as solid lines connecting data points. Further, 
the porosity is printed to the right of each curve of collo- 
cation method results. Results from [Ref. 27] and for a run 
of the Collocation method for variable porosity are plotted 
as diamonds and triangles on Figs. 6 and 7, respectively. 

The results from [Ref. 27] plotted on Fig. 6 are fora 
escullar are airfoil in a choked wind tunnel with no wall 
porosity. The agreement is relatively good but it must be 
considered that the use of a constant fluid acceleration (A) 
is an approximation. Further, the experimental results indi- 
cate a magnitude reduction in the pressure coefficients with 


Introduction of porosity in the manner shown by the results 
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of the computer routine.  Porosities used were 0.0, .05, 
5510.0, and 100.0. 

Figs. 7 and 8 are the results of nonsteady analyses for 
which there is no experimental or theoretical comparison 
"1ο, It can be mentioned, however, that for both the 
increases of porosity produces the expected trend in the 
ΤΙ’ walli tunnel collocation routine results. Note that 
the pressure coefficients shown are normalized to the 
magnitude of the pitch oscillation, αρ. 

It should be noted that for K = .5 (Fig. 8), the Collo- 
cation method solution was near the limit of its applicability, 
as expected from the results of Section III-E. In particular, 


LOr D 0, the source distribution series approximation was 


E 
a Mildly divergent oscillation and the result shown in Fig. 
8-a is oscillatory as well. Interestingly, as porosity was 
increased the character of the results improved. 

mie Collocation routine was modified for the utilization 
ο. πο Following variable porosity, 


σ ον νι το) 


w$ 
iu che porosity varies linearly from zero, at the leading 
edge, to ten at the trailing edge. Of course, there is pre- 
sently no way to verify the results but they are promising. 


However, the source distribution series approximation was 
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more oscillatory for the variable porosity case than that for 
constant porosity, suggesting further analysis is necessary 
Mmemconlirm the Collocation method's applicability to this 


technique. 
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C. DISCUSSION OF MURASAKI AND DRAKE POROUS WALL SONIC WIND 
TUNNEL SOLUTIONS 


Murasaki [Ref. 18] solved the steady sonic flow porous 
wall wind tunnel problem. His results could not provide more 
Bemeone check case against the Collocation method, though. 
However, they show his solution to have good agreement with 
those of experiment for a solid wall tunnel in [Ref. 27] 
(those results are shown in Fig. 6). 

Drake [Ref. 6] Laplace transformed the nonsteady linear- 
ized potential equation and applied the porous wall wind 
ου ο boundary conditions, leading to the following inversion 


for the Green's function, 


gx Tr onen 01.» otprix: sinh[%iXp' 0] | 
3 t e 
ὙΟ1ΚΡ [ ΥΡΚΡ sinhheikp po] +to(ptik) coshfeiky ρ 17) 


πω e iS results of lift coefficient could not be used as 
a verification because of the known inability of tne Collo- 
cation method to produce derivatives of the velocity potential 
for the linearized (zero thickness effect) case. Thus, Drake's 
solution was modified to include the Oswatitsch parameter anc 
computer programmed. The results of that program did not 
agree with reasonable estimates so the Drake solution method | 
Bassfollowes to produce a solution for o IA O EAS 


νά 


test the method at as simple a level as possible. 
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For ge = O e potential is using notation of Section 


III-B 


+ ^ + 
ο = ας) p(X-t.0 ) dt, 
0 


where 


ug = j iK(t - A) 


2 





K^ 
Q o | 
- 40D o] 
O) - $_, / al, ρω 
y-i0o YAP 51ΠΠΝΩΡ p] 
2 
K 
co S X 
= SQ > EEN e n 2 
n=0 pQ 
and 
non 
Sh = a E necem 2. 
p 


προ chis result is derived in a parallel manner to 
Drake's — Z 0, A720 result (inversion by the Theorem of 
Residues), it leads to a series whose convergence is 
@mestionable, i.e., | 


: 2 2 





"m8 
(D 
O 
O 

N) 


29 





in fact, computer generated results using the above for- 
mulation showed good agreement with previous results of Four- 
der and Laplace Transform and Collocation method routines for 
K< .3 and À Z 0. In those cases for which agreement was not ` 
optained, it is believed that the computer routine signaled 
series convergence at a local minimum of term value. 

it was concluded that a Laplace Transform solution of 
this type needs further analysis before it could be used in 


.. τΙΠοπείοη of Collocation results. 
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VI. DISCUSSION AND SUMMARY 


A. ANALYSIS OF SOLUTIONS 
Reference should be made to Section III-E, 
l. Laplace Transform 

Maree olution provided a reliable, if time consuming, 
whe others. Despite the principal shortcoming of 
ata sopaisticated Integration technique, it always gave 
reasonable agreement with the Fourier Transform solution. 
Even for À = 0, the case that required the 'correction inte- 
gral' described in Section III-D, agreement with the Fourier 
solution was good for the parameters tested. Thus, the cal- 
culation of the integrals utilized in the Collocation method 
wuti on was verified. 

Ás programmed, this solution was too extravagant in 
its computation time usage, principally due to the time re- 
“quired for integrating. Presently, integrals are calculated 
from the leading edge to the X value for which a potential is 
desired. Perhaps, a better approach would be to calculate 
the integrals in increments between each X value, as does the 
Collocation routine (see Appendix-C). Further, an improved 
integration routine itself may speed the results and improve 
ΠΕΙ accuracy. 

c xEuourver Transform 
The Hamamoto Fourier Transform method, as extended 


here for variable blade thickness (via the Oswatitsch 
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parameter), frequency (to include K = 0), and interblade 
phase angle, was a reliable and reasonably rapid method for 
velocity potential, pressure coefficient, lift, and moment 
calculations (which includes the solid wall tunnel and free 
jet problems). It must be noted that for large interblade 
Gistances, beyond those required to initially obtain isolated 
wing results, the series involved may not converge. This 
must be considered principally for the A = 0 case. 
bee collocation Method 

The Collocation method, as programmed here, was a 
Qualified success due to its agreement limitations for higher 
frequencies and lower Oswatitsch parameter, as discussed in 
Section III-E. Whenever the other methods showed that the 
potential results were not reasonably smooth or large magni- 
tude differences were present between the real and imaginary 
τ the potential, the coefficients of the series became 
large and the computed velocity potentials deviated from those 
wache Fourier and Laplace Transform solutions. Increasing 
the number of collocation points did not, in general, improve 
the agreement. It is believed that the character of the step- 
wise-linear series approximation in the presence of the numeri- 
cal inaccuracies of the íntegral calculations (though not 
terribly significant in themselves) can lead to large errors 
in the calculated potential, especially near the trailing edge. 

F ο 5 οἱ the inability of the method to reliably 


Podee potentials for à = 0, it cannot calculate pressure 


ps 
O 
N) 





coefficients for that case, which are needed if lift and 
s Ü results are desired. 

Perhaps an improvement to the step-wise-linear series 
is the Fourier series. It might overcome the numerical pit- 
falls encountered by the power series approach of Elder and 
the aforementioned qualification of present series by pro- 
viding a better approximation to the source distribution when 


more complex behavior is expected. 


Ee CASCADE STABILITY 

Two qualifications must be mentioned as to the applica- 
Aent Of these results, produced by the Fourier Transform 
ton. First, the 1 = 0 results were a solution to the 
meari zed potential equation, the applicability of which is 
a function of reduced frequency, K, e.g. Landahl [Ref. 13] 


3 the following requirement for applicability, 


s 


nere, Mr. ΘΓ. ππεπ number. 

Secondly, for t # 0, the approximation of the Oswatitsch 
parameter, A, was utilized. This can only be a rough estimate 
to the real behavior but, at least, it is an initial step. 

A more accurate analysis could perhaps be made through the 


Application of Spreiter's local linearization. (See part 


Momias Section). 
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C. POROUS WALL WIND TUNNEL INTERFERENCE 

The Collocation method applied to the porous walled wind 
tunnel is subject to the same qualifications discussed earlier 
ENEUCHEES Section for the Collocation method. 

The results are difficult to verify experimentally due 
τους paucity of such data. Further, theoretical verifica- 
IDs difficult to obtain as only Murasaki's approach is 
available for steady flow and Drake's analysis for nonsteady 
merous Wall tunnel flows. Unfortunately, Murasaki's data 
could be used to verify only a Single steady flow case and 
the Drake Laplace Transform solution was restricted to the 
imMiamearized case for which the Collocation method solution 
mould not produce reliable results. However, the trends of 
Mco location method solution are correct and the results 
for very large and very small porosities compared well with 
the free jet and solid wall results. 

The modification of the routine for variation of porosity 
along the wall appears promising and may allow a numerical 
Pest ot the Sears [Ref. 25] technique. Further, a kind of 
MOCcal linearization may be applied to collocation methods 


Meee tne next part of this section). 


D. FUTURE WORK 
Various areas of improvement have been recognized. 
1. Improve series approximation methods for Collocation 


methods. 


104 





2. Improved numerical integration schemes for those methods. 
3.) Abplication of Spreiter's 'local linearization' to 
Fourier Transform and Collocation methods. 

The first two have been discussed. The third requires 
της following elaboration and description for steady sonic 
flow. 

Recall that the Oswatitsch parameter was used to approxi- 


mate the nonconstant coefficient of the governing equation, 


νε. 
ὀγὺ - (ΥΠ) ὀχχόχ = 0, 
Became, 
yy ~ Ardy = 0; 
where, A = constant > 0, assuming that the fluid acceleration 


was slowly varying. 


mess olution for the steady flow potential on a single 
KOIL is, 


naso Ne nr 


VTA =) ER 





Now, Spreiter argued [Ref. 25] A can be replaced by its 
definition, yielding an ordinary differential equation whicn 


admits an analytical solution. 


TOS 





n τος recently extended this approach to the oscillating 
@erroil in transonic flow (Ref. 28]. Start with the oscillating 
INI Potential equation for sonie flow, 

2 
ES) 
Hi i α 5 
Y(X,0) = - —— { v(s) ———— QS, 
WA -1 X-S 


where, 
Q = À + 21K, 


À = constant, replacing the variable coefficient (+) dy 
With y(X,0), calculate V4 5,0), replace A by the function 
it Originally represented, and solve the resulting second 
order differential equation for y(X,0) and v„(X,0). 
It remains to be investigated whether this approach could 
be further extended to study nonsteady interference in trans- 
Ms. Tt appears from the outset, however, that the 


Eouocatdion method is suited to a form of 'local linearization!'. 
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Vit.  CONCDOSTIOSIS 


The unsteady, sonic, small disturbance potential equation 
was approximated, following the approach of Teipel [Ref. 30] 
and Hosokawa [Ref. 11], by replacing the nonconstant coeffi- 
cients with a constant, the Oswatitsch parameter. The re- 
Sulting parabolic equation was solved by Laplace and Fourier 
Buusrtorm and Collocation techniques. 

Comparison of solutions to numerous sonic cascade problems 
usted that the Transform solutions were reliable and con- 
Lent in the Fourier en Case. ine  Cellocary tons vech— 
nique was usable with qualifications, growing out of the series 
approximation and numerical integration methods utilized, 
which degraded its accuracy for certain parameters. 

Son+e cascade stability boundaries were computed utilizing 
the Fourier Transform solution and interference was found to 
be strongly destabilizing. Such results must be qualified by 
the realization that they were obtained by the solution of an 
meproximated equation, though the trends are not expected to 
be materially affected. 

The Collocation method, once some confidence was obtained 
from its comparison with the other solution as mentioned 
above, was applied to the sonic porous wall wind tunnel inter- 
ference problem. Though few experimental results were avail- 


able for comparison and the Drake [Ref. 6] approach was 


LO 





ου ioned, the trends of the results for different values 

of porosity were as expected and solid wall wind tunnel and 
sonic free jet results were easily obtainable with reasonable 
porosity values, 

Miiemsuccess Of the Collocation solution in handling a 
range of sonic interference problems leads to the considera- 
tion of improvements beyond better series approximations and 
numerical techniques. It is reasonable to expect application 
Of the Spreiter [Ref. 28] ‘local linearization! to improve 
accuracy. Further, the method appears to be capable of uti- 
lizing a streamwise variable porosity in a manner like that 


proposed by Sears [Ref. 25]. 
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APPENDIX A 


RORMUEATTION OF THE LAPLACE SOLUTION ROUTINE 


A-1 Main Program 
The Main routine serves as the controller for calculation 


of potential utilizing the Laplace Transform solution derived 


in Section III-A. x 
K* 2 
ο. τ. | = (2nptp) 
b (x ,0") m J mn... αὐ; (vf S) e"? [2e (X78) 4 
47Q -1 ES  n=0 
_Q (2np+2p)“ Q (2np)* 
- vo (S)[e UPC uu O SS. πι ο) 


Calculations are made ο eontr ronne ya tables DK. ΠΠ 
RHO, and SIGMA (corresponding to K, à, p, and o in equation 
ΘΠ) at each of L selected points along the chord. 

To reduce computational redundancy, the Main routine 
calculates the parameters A, ©, and D as described in Section 


III-D. Also, the boundary condition expressions 
vg(S)e ^^ , 
and 


va (5) 
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are treated as they are in the Collocation routine (Appendix 


C) 


The Main routine reduces equation (III-22) to 


where 


and 


v (s) 


1O 
e vo (S) 


POT(J) 





O 
-l - iKX 
| 


R * I, 


2 


η. 


PIO 


(A-1-1) 


(A-1-2) 


(A-I-3) 


(A-I-4) 





Q (2np*p)* 
T 


2 
oo X Exp[t-(X-8) = πα ο 
m 2 {5.6.0 A = 
mu ΕΞ uYX-S 
2 2 
X Exp[-—(X-8) - i des 
+ 23611 το 
ET 4x-S 
2 
Kau Q (2np+2p)° 
SPP O (YO, ys 1 
E: 90 ME 65 
E ES 
9 2 
K Q (2np) 
X Expl —(X-5) - 
Jg: f mu hon ERES X) ως ] 
E X-S 
i X ert ia τ Q ESL on 
. zn [f S Bc M τ. 
E NS 
X τω ees „a — 
"ES o. 9 È EST! 481), EOS 
E JS . 


which is calculated by subroutine ISUM. 

The Main routine uses the common named SUM to relay 
boundary condition and integral sum results between itself 
Emdesubroutine TSUM, ¿Common HOLD relays parameter information 


between Main and subroutines ISUM, ZSUM, and INTCOR. 
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A-II  Subroutine ISUM 

ISUM is called by Main for each X position on the blade 
Teh the potential is to be evaluated. It sequentially 
113565 each term in the summation comprising I in equation 
(A-I-5), in the 'n loop'. 

ISUM denotes the integrals in equation (A-I-5) as the 
s les Iji where j = 1,2 and i = 0,1,2. Since ZSUM is 
ucusror the evaluation of the integrals, jJ and a constant 
denote the power of U required in the integrand (note that 


this routine's ZSUM calculates the integrals with uJ”1'> 


> 


vice y4ré. Si us treated in section Ill-D). 
Because each integral differs not only in j, but in the 
value of the natural exponent values, the exponent is 


considered to be of the form 


τε = 2np, Rl = 2nptp, R2 = 2npt2p. Thus, once in the 
'n loop'!, ISUM ealeulates RO, R1, and R2 and the associated 
Bi and Ei required for the subsequent integrations. 

The integrations are conducted sequentially on i. Since 
NP 3nd IJO for n do not differ, integration for IjO 
can be avoided for n > 0. 

In the event that à = 0, correction integrals must be 
calculated as outlined in Section III-D. Subroutine INTCOR 


is called for the calculations which are stored matrix 
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COR(i,jJ,k). The corrections are to be calculated but once 
Borzeachn; thus, NMAX is set on the first call of ISUM and 
subsequent calls use it to determine whether further correc- 
tion calculations are required (resetting NMAX higher). 
Accuracy of the correction integrals is dependent upon 
the magnitude of the ratio (E/U)/DU, where U takes on the 
ο value used in the correction integral, i.e., U=XA 
Meeemocction III-D). If we desire (E/U)/DU > 10,000, since 


Bea k= 0.1,2, 


MLA 


then, the maximum value XA can be 


XA = (2n t k)p/ 100. 


O AS calculated for each n and i. If the calculated 
Ae Ss oben the required X0, it is used. If it is greater, 
an XA iS chosen smaller than XO. 

p s rr e for 4 7 0 and RO = 0.0, XA is chosen as 
1.0D-03 and ZSUM calculates a correction (see Section III-D). 
Subsequent integrations utilize the arbitrarily small 
XA = 1.0D-08 as a compromise between computation time and 
accuracy. 

Convergence of the series for I is signaled wnen successive 
terms alter it by less than DEL x 100 percent and control is 


returned to Main with the value for I. 


jll. 





SUBPOUT INE ISUM(X) 


LERAKZOO) 


DIMENSTON 


O 


O 


ας ο 
ON e 
Ç) O 

== ji 
DIN oem 
n" 
bh! e” 
EU ru 
LEE fr 
uc uu) 
NON 


Qt, 


XO 


`. 


.ODO 


DO? 


e 
O Owi 
O + Oe 
= ex 

οὓς li H 
Ol HOI 
iu Ou 
e XXX 


Cy CONJ 
€ C. CY O — (J 
st 3° cx Y x 
Ormnss 3: 3% 
LADON 
3⁄8 3€ 3° CY C: _ Ωα 
ZA O SF 
— T T C. O O: >x >< `< 
F- c uc OOO 
«I O O C +z + t 
C^ 4- + CO Coco O 
ri J WLVLA WV OO 
| LL OHAINNINMIWN 
OR Cr ë 9999 
IT AT 
Co lI QCN GO CN GO =+ CN 
CZ vr e eco cou αἱ 


LOOP. 
N1=1,L IM 
)*RHI*2.0D0 


109 


e 
< 


THE N 


INEO 


«ii 


OO 


am pr 


CO 

"(N 

Im 

-~ o 

AAN] 

~ nl) 
OON 
XX 

° e 

«I« | 

>< >< 

er er CO 

= Z.— 
am 
UI MN 
NNX m 
— gg dd 
AIO O 
{κους ντ 
OORO 


N 





120=122 


+ 

O 

EL 

O 

c^ 
° 
e~ = 
E) vU» 
e a. 

LLJ 
Η! 0 
= ΕΗ 
<I C5 
ΚΠ 5 
= 26) 
ο τι) 
Du 
<í H H UL 
C) GO LU e 

e 
VO 


O 

O 

ο 

0 

O 

O 

ez 

NS 

C: LA 

m 

α O om 

eI Ea 

— >< >< 

Wu O e 
A OY 

+ <I e 
x ur 

Z >x Z 
coe < we 

NC <a 

=> ZO 

— ° e C 

«ful W- 

EJ wa Y 

el Crd Ors 

Uu «ισα (0) 

OXXZ Jr 


od ft mm” |} me 
uU < LL. <1 U. <í G 
NS >< — >XX OD 


e e e e 
HANAN 
PP HA ERA 
ο» e 
του + + 
e e A 
Οσο) C 
A AAA AA 
94 vr, + e 
«(fcc | 
>x >x Z Z 
CO wer ο *.— 
OSs >= HANIA 
Quoc o + 
ONNC OX 
° LO 
rA dJ- JH HH 
ved edt — 
«1 «1 «1 —OC 
AX CQ —— 


TU 


D 
(9 
= 
O 
co 
e 
CN ane 
RP Gi 
p Q. 
LL: 
UJ 6 
= I- 
<I C5 
JJ 9 
=æ) =] 
ONNO 
-s Œ) U i 
«τι. 
O Q US a 
C CO 


ή 


c 
O 
e 
e 
e 
O 
e 
N 
C «ο 
s p 
c Om 
HA [-- 
mx >< 
Wu Cie 
eI ον 
+ «A e 
>x e 
= >Ç Z 
J} -- < Í w 
NO > m 
Ἂς ZO 
fp ο ec) 
Iu U 
ETS a 
wt LO] te 
La ar 
Ox <2. 
Nittel J 
|! «iu. «ilL «i 


x >< — >< —: COD 


— om 
ο ο οι ον 
ec C 0C 
bad fd Od I 
” ° 
eN + + 
ο ο Oy 
OO meV 
> XX NA 
e e... o 
<] AAA | 
>x > Z Z 
wrr p eN) 
ΞΘ «Γι 
Dune: 
t u yC LK 
-ω-αζς ΟΚ 
De) 
O α 
d e NENN 
< <1 — CN (NJ 
COC) bt bt tee 


ο 


Ow) 


TERT 


OC 


]J22+120)) 


(THOO*(T112+110) + THO1*( 


ASNO 2} 


OD 





N3 

Νας 3 

JüNILINOO 

Ναι 3 Ὁ 

QOI οι 09 {Π30519”λΗ93}31 
(1/land3i - I))SE8VQ9=HHD 
20-00*T="136 

*NOIY3i139 1LN329293d *T 3Η 35η 


001 


o 


l" 





BENT Subroutine INTCOR 
ames OY TSUM in the case A = 0 INTCOR follows the 
organization of Elder [Ref. 7] in calling the corrections 


required for the integrals from the subroutine COR. 


A-IV. Subroutines COR and ZSUM 

eons a Simplification of the more general correction 
routine developed for the Porous Tunnel Wall Solution. ZSUM 
is an adaptation of Elder's [Ref. 7] pyramided Simpson's 
Rule integration routine with P = 0 corrections added (see 


See ion III-D). 


A-V. Subroutine CS 
This is a orar routine neil Izcd to provide Fresnel 
Integral values required in À = 0 correction integral 


calculations. 
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SUBROUTINE ZSUM( XA,X?2 4, J4 SUM) 
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APPENDIX B 


FORMULATION OF THE FOURIER SOLUTION ROUTINE 


B-I Main Program 

The Main routine satisfies equation (III-47) derived in 
"enon T II-B for potential at Selected locations, X, along 
the unit blade chord and equations (III-49, 51, and 55) for 


pressure coefficient, total lift and moment coefficients. 


y(X,0) = -((C-1-1K, A.) [D. T. (X) €B: 8, (X) ]-1K, LD- T5 GO *B: 55 00 1] ; 


(EII-47) 


Cp(X) = -2[ y, (X,0) 5 ik, p(X,0) ] : πο 


e 
I! 


Z9) s iK F4(1.0)] ; (21157) 


C3 
l! 


-2(1K, [F,(1,0)-P3(1,0) J+[9(1,0)-F,(1,0) ])-A, *C, > 


μυ 55) 


where Tw and ὋΝ εν" viqeqhoy P up rontinmerSUM as describes 

che next section and En are successive integrations of y. 
Main controls the calculation of TN and DN With Subroutine 

SUM in peculiar manner. Since ὋΝ and TN arelaecnually iniinite 

series, Main calls SUM to produce successive terms in those 

ο ο ο the sums of which are ὋΝ and Tu: 


Convergence is signaled for the loop of index n in Main 


when successive calculations of the potential, for a Dare euler 
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Boerse posit*on, differ by less than criterion EPS. 

When the lift and moment are calculated, the n loop is 
memmrnued, if necessary, to provide the successive calcula- 
tions of the magnitude and phase angle of the moment coeffi- 
cient until convergence is detected in them. 

Due to the difference in nondimensionalizations between 
Pee Fourier and other solutions, RHO, DK, and DL must be 
Beleg, as shown in the output, to conform to the solution. 
IM ner, the results for the potential obtained will be one 


Raul Ghose of the other solutions. 
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ROUTINE USING FOURTER TRANSFIRM APPROACH TO CASCADE 
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Beil, subroutine SUM 
SUM is called by Main, as explained, to calculate the 


memes for S. and T... Recall that S. (X) and T- (X) are from 


N N 
equations (III-66 and 67). Main provides the coefficients 


to the series as Z. 


Since 


X 


E So) dg 
Syed = B N 


and likewise for Tw» then the dummy variables Un in SUM 


provide the appropriate series terms for Sy or TN for the 
value of J (J = 2 for Ty? sent to it in the calling statement. 
Meses provided to properly caleulate the first term 


of the series for T. for the case K = O, 
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APPENDIX C 


FORMULATION OF THE COLLOCATION SOLUTION ROUTINE 


C-I. Main Program 

Because the derivations have certain parallels in the 
integrals to be evaluated, much of this routine and the 
Laplace Transform routine are similar. 

Main calculates the necessary variables required by the 
integrating routine, the boundary condition for coefficient 
and potential generation, and the step size between colloca- 
Bon points (XINT). With that done, Subroutine INTEGR is 
called to calculate ali integrals required in subsequent 
Sgeps. Main obtains the coefficients from Subroutine COEFF, 
then, implements equation (111-85) and calculates the potential. 

Men's organization of the solution for the potential is 
Also applied in COEFF and is due to the series chosen. The 
J loop carries the calculations down the chordwise collocation 
points and the I loop conducts only those calculations as 
allowed by the chord position dependent step function in equa- 
Bone 1-5). For any X position, at least two terms must 
η Ἰατεὰ, the n= 0, 1 terms of the series, and are done 
when I = 1. For I = 2, no calculations need be made and it 


iessidesteppved and I is cycled to 3. 


1 





AND SIGMA. 


DK,DL,RHO, 


06060 00509000000000500000000000000000000000000000000000090000000000000.. 


INITIALIZATION SECTION 
SET THE CONTROLLING VARIABLES 


COMMO 


COMM 
CCMMO 


ως) ως.) 


σι ο 
"COO 
ADO 
NON ο 
TÍ] a ei 
> O O !! 
On wo 
— x lI 
ως Qua 


Φ 9 @ e 9 9 Φ Φ @ 9Φ 964 Φ @ @ @ @ @ @ 0 @ @ @ @ @ @° @ @@ @ @ Φ 9 9 @@ @ @ 499 


O 


i 


^ 


4 


Gi 


@ t--t 
ο.) 


0 y” 


. e 


9.0 
@ ur 
eu 
εἰ-- 
ere 
e° r 


o > 


-— Am 


À OK OL, RA 


ABLES THAT ARE 
ON VARIQUS 


V 
DD 


FIXED VARI 
IN ADDITI 


a ° 


OLD GOWOU 


N 

3% 

>£ 

Y 

CA 

3% 

O 
ra A 

OO O 
C eC) = 
|-- + O 

+ 0 e 
C3 cN ev 
(5 e ~ 

AN AN 
~ E x3 
O OD + m SC 
t ~~ C) OO 
> ~ = = 


SOONG (Ey 
e> + * {Ot 
MOD eas 
O sO ‘DO 
O Mm ON 

u- HCO H HW H Ú H 
C «1 5 <( CQ (O CO) LL 


aN 


DCMPLXCOL,2.0DO0*DK) 


e 
== 


Q 





(N*€)1NI — OlIdxX =IIId 
(W*Z)INI -Olld 
jOv18 iN39VlüV 

(W*S)1NI — OOIdx*X =TOId 

(W YILINI -OO0ld 

Jüv18 3oN3x3J3v 

I4 I-P-W 

+ 0l 09 ο ο. 

C*t=I > DO 

090 * I-=NX 

πέτ-ρ ; 04 

(0:14) 140S0I/000* T-=3 

000*2/1NIX + 000*T==X 

*SIVIINILOA 3HL YOY 3A109S 
(Ns τα; 
/*,0152)0N) *X22* 1(I* DON. *XS*.. ΞΙΝ3Ι31344303 SHI .*/7/)ivWuod 
(1351s I* e I*2)0NSLI* T)AIN). (9006* 9) 31i1uM 
(IŚiNIX)33309 11V9 


το ΝΗ “ΙΗ: O NOA ONO 


%6..0.00000.0.0un....0.1n.0.0..EOe. u, .„“s—...eae„aaaeeaua„..,..,a„aeus„esear sr ...s„....e„.ea eu ee 0 so 910 0 8 8 z os 2 oo en oo 


(1* 12NIX)9931NI IV) 
"NOL1AVINITWO οσο” SNA Tuy MV 


99499 99 Φ Φ Φ Φ 9Φ9Φ 6 99Φ 66 Φ 99899 ΦϑΦ 98ο 99 6898994996 96 9 6ΦΦΦ9Φ6Φ 96969896 ϑσ“ 996 8Φ6Φ6Φ6Φ599949Φ99Φ 6Φ6Φ 69599 9 Φ 


(1)1v0130/0600* 


IVAYJLNI ONY ‘JNO NOILVLS “7103 3H1 NO SNOIAVIS 30 33I8WN; 
TOH1:9 I 

OOH 9 I 

(VWoOISx(OGO*T1*000*0)X' IaWw2a)dX 

(0000-2 en ΠΤ 

"SNOILIONOS AYVONNOS JHL 3.170 

H 


1 
0 
I 
I 
0 
I 
(14/0)1305039x*(000* 2 /OHY 


QOO 


zi 
WOW 


7006 


OVVVO OOIOOODOO 


QUO 


OW 


195 





( (OTIdxNX 


TT1d)x(I° 2) AN 


1721025 


(OOId*NX 


1041. t ° 2 TG: 


"uz ια zoo 


za οο νο = 
xD Zu.xonNUl 


- TOId)x(I*T)ON)«u * 10d -10d 
"SWu3l ININIVWIY JHL MON 


UE) 
LNIXXOGO*C + NX=NX 
THL) = OTIdEl(I ZIAN rot 
00Id*((T*TINN + 00H41) )xuziOd 
*“SAY3IL LNVISNOD 5Η} 15814 
SEO ZI Et 


+ ON ~ 


Or 


OVO O 


OVO 


136 





C-II.  Subroutine COEFF 

This subroutine is called by Main after all integrals 
have been calculated for the purpose of solving equations 
(III-86 and 87) for the coefficients of the series representing 
Premenverference source distribution. 

For n = 0 and I, the abovementioned equations can be 


I ameed to matrix form, 


. - 1 X RIO RII WR -BBO 
= E 
"πο ΠΤ il X vo] I -Bnl 
1 X RIO RI1 το -BBO 
n=1 [ 
RIO RII il x J I |-BB1| 


ine first two rows are calculated at the first collocation 
point, selected as X = -1.0 + XINT/2.0, and the next two at 
po ccond point, X = -1.0 t 1.5*XINT. 

The above matrix equation is solved by a compiex matrix 
inversion routine (CMINV) and a multiplication to yield the 
Hscusolents NU(I,J), I,J s 1,2. (1.e., Voo2%o12Y192 2nd νη): 

The presence of the step funetion in the equations greatly 
simplifies subsequent calculations. A matrix solution is 
used again but only a single left side matrix need be calcu- 
lated and inverted. As successive collocation points are 
ου ο οἱ previous coefficients have been calculated and 


they are used to calculate a new right hand side vector. 
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CALL CMINV(AA+,2+D+LL,MM) 
MULTIPLY FOR TAE RESULT 


RETURN 


SOLVE. 
END 
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1. Subroutine INTEGR 

eause OF the variety of integrals required, they are 
calculated as soon as possible and stored for further refer- 
ence in the matrix INT(5,m) carried by Common INTG. 

Since the integrals are of the same form as those encoun- 
tered in the Laplace solution, the same solution is used, with 
Manor differences for the differing exponent J range. Sub- 
routine ZSUM, borrowed from Elder [Ref. 7 ] and modified for 
ου 0.0, is utilized for the majority of the integrations. 
In the event that the Oswatitsch parameter is taken as zero, 
ο cine COR is utilized to calculate the ‘correction 
mauecralt , 

MER provides five different integrals, for the following 


ο ο οἳ P and J values, 





J 
L-1 
for the ranges 0 to X = -1,0 + XINT/2.0 + 2 XINT because 
n=2 
X X A 
ΝΤ. 5/45 = ο ο des. f íf (y 
0 τ X 0 
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where 


ο οσο, the integrals are calculated stepwise, i.e. the 


range of integration starts as O to X ioc OM οσο, 


ο. 0 1? 
until all necessary integrals are obtained. 
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SUBROUTINE INTEGR(XINT,L) 
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IF(B.GT.EPS.-OR.P.EQ.0) GO TO 3 
DL=0 
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C-IV. Subroutines COR and ZSUM 

Subroutine COR is an adaptation of the more general 
ΙΕ τον the 'correctíon integral' used in the Porous 
Tunnel Wall program described in Appendix D-V. 

Subroutine ZSUM is a modification of Elder's [Ref. 7] 
pyramided Simpson's Rule integration routine to include 


p = 0 corrections (see Section III-D). 


| 5ubroutine CMINV 
Subroutine CMINV inverts the complex matrices for 


Eusutine COEFF. It is an outgrowth of a Library Subroutine. 
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SUBROUTINE ZSUM(XA+XOD+ Jy SUM) 
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APPENDIX D 


FORMULATION OF THE SONIC POROUS WALL WIND TUNNEL ROUTINE 


D-I. Main Program 

This Main program is very similar to the Main program 
1η Appendix C-I for the Collocation solution of the 
ου problem. Provision is made in this routine, however, 
for the placement of the 'tunnel wall' at Y = p/2 (which will 
be RHO in this program). Furthermore, two different Oswatitsch 
Pasa ceps may be utilized, one for the airfoil, DL, and one 
for the wall, DLW. Because the wall has no ‘isolated blade' 
ο. distribution, THIO and THll do not appear. Also, đe- 
rivatives are made available so as to calculate the pressure 
pouwesdent as well as the velocity potential. Finally, the 
variable, SIGWL is the porosity factor (see Section V), which 
can be made variable along the wall (see Subroutine COEFF 
here). 

Meson program not only controls the calculation of 
Ma teerals required. by calling Subroutine INTEGR, and 
obtains the appropriate interference source distribution 
series coefficients from COEFF, it calculates the potential, 
ο -υτ ας the pressure coeíficient, CF, at the selected 
ios ons (collocation points) along the chord. 

From Section V the potential on the airfoil is 


P(X,0) = $9(X,0) + 0,(X,0) + 6, (X,Y 2-0). (D-1-1) 


153 





Also 


O de 
Y, (Xx,0) EE Pox (KO) T oy (X50) + Py Y =P). (D-I-2) 
The pressure coefficient is 
Cp(X,0) = -2[Y,(X,0) + 1Ky(X,0)]. (D-I-3) 


By introducing the step-wise-linear series used in the 
Ful ion solution for the cascade (equations III-81 and 


82) equation V-2 becomes, 


L-1 
= + 3: n iE 3 
y(x,0) Ro {Lepo 99g 1P100 + [θρ1 νρι 4101 ο. ΟΝ } 
L-1 
4 1 mm 
*R4Lv, 9PI10*v, ,PIIl Coin E , (D-I-4) 
where, 
X X X 
P100 - / Ees, prov =x 2L ,HEXÑplcl_ as, 
-1 -S S ES =] (X-5) 
Ke 
PION = PIO1 - XN-PIO0O, C = a SE 
0 


R2 - p/2- ο (al i) 
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and, for PI10, PI11, and PIIN, 


2 Q, 2 
ὦ“; 
A 


ο ος e: 


3 
K or The X derivative of the potential, all of the 

Xx terms of the potential are preceded by a D, denoting 

the derivative. The appropriate differentiated integrals 


ape obtained from Section V-A-3. 


Ke Ke 
Exp{—(X-S) ] Exp[——(X-S )] 
3 x Qo Qo n 
SS eS 
E Nm qx 


etc. 
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Or. subroutine COEFF 

As the Main routine, this routine differs only in certain 
Bl su from the Collocation routine for the cascade. Of 
greatest importance is the handling of the wall boundary 


condition equation, equation (V-5) 
doy (X,Y 70) + doy (X,Y =P) + 9, y(X,Y,=0) + 
0 c | = z 
O, 9 Hoy (Xx Yo07?) + Poy (A Y 70) t $44 05, Y,70)] 
* iK[$9(X,Yg70) * 6g, Yo) * $4 051,70)1) 7 0. (D-I-7) 


NU Corsanization of the solution for the coefficients is 
ΕΠ “5 that in Appendix C-II save for the inclusion of 
the SIGWL, the deletion of the isolated adjacent blade poten- 
(Ian normal velocity component, and a slight change in 
M@mation. 

meimstrated is the approach for varying the porosity. 
SNS calculated by the appropriate function in COEFF and 


used as necessary. 
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ITT subroutine INTEGR 
This routine is similar to that described in Appendix 
C-III except it is more complex to calculate the additional 


integrals required. 
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E Subroutine COR 

Subroutine COR calculates the 'correction integral! 
when commanded by Subroutine INTEGR, using equation (III-99). 
ities correction integral’ is obtained by summing the terms 
Dor uncated series, here SI(I) and CI(1), I = 0,1,2, and 3. 


SI(I) and CI(I) are integrals of the form 


ΠΡ ym/2 cos 


: G [I 
which can be evaluated by integration by parts and use of the 
Fresnel Integrals (see Section III-D). Elder [Ref. 7] noted 
that this is possible and reasonably accurate for powers of U 
Weseeunan about 1.5. In any event, integrals of larger powers 
are very small. 
COR calculates the terms by liberal use of equations 

(III-101 and 102). It must be noted that for j = 0 the use of 
sn or (111-101) is not direct, i.e., j = 1 is inserted into 


(111-101) and the resulting integral on the right hand side 


is the desired j = O integral 
/ P “AIBN sin παν G 
0 


Thus, it can be seen why the involved logic is contained in 
COR so as to calculate ocrrections for j = 0 without destroying 


miembasite form of COR. 
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G course, once j = 0 and 1 corrections have been 
calculated no further calculations of terms must be performed 


as they are deemed adequately small to neglect. 
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AND COS TRUNCATEC SERIES TO OBTAIN THE FINAL 


HE SIN 
CORRECT ION. 


> 


OOo 


w wa SS 
pd or (C 


uu oe 
αο u D 
Cca 
OO OW 
MO Nu) 





DEE »ubroutines Z5UM, CMINV, and CS 

mese suorolutines are the same as those used in the 
Kise Lion routine for the cascade. The ZSUM routine 
Fr n ποπενες. in that the terms required for the 


'correction' in the case P = 0 are precisely equations 


eens), to provide greater flexibility. 
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